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APRESENTAÇÃO 


O presente texto traz os resultados mais importantes obtidos a partir 
do trabalho de pesquisa desenvolvido pelo autor, durante a elaboração da sua 
tese de doutorado, junto ao Centro de Estudos de Engenharia Civil - CESEC, 
da Universidade Federal do Paraná, junto ao programa de Pós-Graduação em 
Métodos Numéricos em Engenharia. 

O tema principal desta obra versa sobre métodos de lagrangeano 
aumentado, os quais são conhecidos, populares e estudados pela comunidade 
acadêmica que trabalha na área de análise numérica e otimização. Boa parte 
desta popularidade do método, se deve a sua grande aplicabilidade na resolução 
de problemas práticos ou teóricos, principalmente nas áreas de engenharia e 
na matemática aplicada. Tais métodos têm em sua estrutura, como componente 
principal, uma função (termo) de penalização. As diferentes variações dessa 
função de penalização podem dar origem a diferentes métodos de lagrangeano 
aumentado. Seguindo esta metodologia, neste trabalho apresenta-se um 
novo método de lagrangeano aumento cuja função de penalização é do tipo 
exponencial. Apresenta-se o modelo, algoritmos relacionados, resultados sobre 
condições de convergência e implementações. Por fim, esse método é aplicado 
na resolução de uma classe específica de problemas matemáticos importantes: 
problemas de equilíbrio. Nesse contexto são apresentados resultados numéricos 
referentes às implementações dos algoritmos, análises, conclusões sobre os 
resultados e, também, apontados caminhos para a continuidade da pesquisa. 
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RESUMO 


Neste estudo propomos um Método de Lagrangeano Aumentado Exponencial para 
resolução do Problema de Equilíbrio Geral. Tal método é uma extensão dos algoritmos 
apresentados em (NASRI, 2010). Em nossa proposta, substituímos a penalidade clássica 
de Rockafellar presente nestes algoritmos pela penalidade exponencial. Em seguida, 
refizemos a teoria geral em torno do novo algoritmo. Mostramos que o mesmo está bem 
definido, que possui terminação finita, que admitidas às condições iniciais adequadas 
o mesmo gera uma sequência convergente. A teoria de convergência é construída via 
ponto proximal, ou seja, a partir da relação de dualidade entre lagrangeano aumentado e 
ponto proximal, porém o termo de regularização aqui utilizado é uma quase distância de 
Bregman, diferente do termo quadrático empregado em (NASRI, 2010). Para contornar 
possíveis problemas de mau condicionamento causados pela presença da penalidade 
exponencial, procedemos a um ajuste quadrático da mesma. Em seguida, testamos 
a nova metodologia por meio de experimentos numéricos, considerando inicialmente o 
método com a penalidade exponencial e, em seguida, com a quadrática ajustada. Para 
realizar os testes, escolhemos problemas de equilíbrio e problemas de equilíbrio de 
Nash generalizados (GNEPs), os quais compõem uma classe particular de problemas 
de equilíbrio. O método puro resolveu problemas pequenos e, para problemas com 
dimensões maiores, a versão com quadrática ajustada mostrou-se melhor. Na sequência, 
discutimos os resultados numéricos e apresentamos as considerações finais. Para 
finalizar, deixamos algumas perspectivas de trabalhos futuros e uma lista de referências, 
as quais serviram de suporte para essa pesquisa. 

PALAVRAS-CHAVES: Palavras-chaves: Problema de Equilíbrio, Lagrangeano Aumentado, 
Penalidade Exponencial. 


Resumo 


ABSTRACT 


In this study we propose an Exponential Augmented Lagrangian Method to the resolution of 
the General Equilibrium Problem. Such a method is an extension of the algorithms presented 
in (NASRI, 2010). In our proposal, we replaced the classic penalty Rockafellar present in 
these algorithms by the exponential penalty. Then, we redid the general theory surrounding 
the new algorithm. We show that it is well defined, has finite termination and, admitted the 
appropriate initial conditions, it generates a convergent sequence. The convergence theory is 
built via proximal point, ie, from the dual relationship between Augmented Lagrangian and the 
proximal point, but the regularization term used here is almost distance Bregman, different ftom 
the quadratic term used in (NASRI, 2010). To work around possible bad conditioning problems 
caused by the presence of the exponential penalty, we proceeded to a quadratic adjustment of 
it. Next, we tested the new method by numerical experiments, considering initially the method 
with the exponential penalty and then with the quadratic adjusted. To perform the tests, we 
chose equilibrium problems and Generalized Nash equilibrium problems (GNEPs), which 
compose a particular class of equilibrium problems. The pure method solved little problems 
and for problems with larger dimensions, the version with the quadratic adjusted proved to 
be better. Following, we discuss the numerical results and present the final considerations. 
Finally, we leave some perspectives of future works and a list of references, which served as 
support to this research. 

KEY-WORDS: Equilibrium Problems, Augmented Lagrangian, Exponential Penalty. 
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INTRODUÇÃO 


11 MOTIVACÃO 


Em anos recentes o problema de equilíbrio tem recebido atenção dos 
pesquisadores. Diversos resultados teóricos e experimentos numéricos têm sido 
publicados. Entre os trabalhos desta natureza podemos apontar (BIANCHI; PINI, 
2005), (BIGl et al., 2013) e (BLUM; OETTLI, 1994). Ainda, sobre esse tema destacamos 
(BURACHIK, 2012), (FACCHINEI, 2003), (IUSEM; NASRI, 2007), ([USEM; KASSAY; SOSA, 
2009), (IUSEM; SOSA, 2003a), (IUSEM; SOSA, 2003b), (MASTROENI, 2000; MATIOLI; 
SOSA; YUAN, 2012; MUU; OETTLI, 1992; MUU; QUOC, 2009; FACCHINEI; KANZOW, 
2007; QUOC; MUU, 2004; TRAN; DUNG; NGUYEN, 2008; SANTOS; SCHEIMBERG, 
2011; NGUYEN; STRODIOT; NGUYEN, 2009; NGUYEN; STRODIOT et al., 2009). Em 
particular (BIGI et al., 2013) apresenta um apanhado sobre algoritmos e condições de 
existência de soluções para esta classe de problemas. 

Como casos particulares do problema de equilíbrio destacamos: problemas de otimi- 
zação convexa, problemas de ponto fixo, problemas de complementaridade, problemas de 
equilíbrio de Nash, problemas de desigualdade variacional e problemas de minimização 
vetorial, como destacado nos trabalhos de (BIGI et al., 2013; BLUM; OETTLI, 1994; IUSEM; 
SOSA, 2003a). 

Em (NASRI, 2010) são apresentados métodos de lagrangeano aumentado 
para resolução de problemas de equilíbrio em dimensão finita, cujos conjuntos viáveis 
são definidos por inequações convexas, generalizando o método de lagrangeano 
aumentado proximal para otimização restrita. A cada iteração as variáveis primais são 
atualizadas resolvendo-se um problema de equilíbrio irrestrito e, então, as variáveis duais 
são atualizadas por uma fórmula fechada e a análise de convergência é apresentada. 

Na sequência definimos o problema de equilíbrio com as propriedades mínimas 
necessárias para o desenvolvimento teórico. 


Definição 1.1 Considere K um subconjunto de Rº não vazio, convexo e fechado 
e uma função f: K x K — R satisfazendo 


Pi: f(x,x) =0 para todo xe K, 
P2: f(x, ):K> Ré convexa e semicontínua inferior para todo xe K, 
P3:f(,y):K> Ré semicontínua superior para todo ye K. 


O problema de equilíbrio, denotado por EP (f, K), consiste em determinar x* e K 
tal que f (x*, y) > 0 para todo y e K. O conjunto de soluções para EP (f, K) será denotado 
por S(f, K). 

Em termos de métodos computacionais para problemas de equilíbrio destacamos 
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os trabalhos de (BIGI et al., 2013; FACCHINEI; KANZOW, 2010a). Em (NASRI, 2010), 
também são introduzidas versões exatas e inexatas de métodos de lagrangeano 
aumentado para resolução do EP (f, K) em Rê, para o caso nos quais o conjunto viável 
Ké da forma: 


K=xeR :ho)<0 (I<iem) (1.1) 
em que h,: Rº> Resão funções convexas para todo i= 1,...,m. 


Observe que na Definição 1.1, há a necessidade das duas variáveis, x e y, 
estar no conjunto convexo K, ou seja, resolver o EP (f, K) consiste em determinar x* 
ekKtal que f (x*, y) = 0, para todo y em K. 

Devido ao fato das variáveis pertencerem ao conjunto convexo K, a maioria dos 
artigos encontrados na literatura recorrem a algum tipo de projeção para atender essa 
condição. É bem conhecido que fazer projeção resulta em outro problema de otimização 
que pode ser um problema de difícil resolução, a menos que o conjunto K seja tipo 
poliedros. 

No caso dos métodos de lagrangeano aumentado, a “dependência de estar no 
conjunto” será enfraquecida através da abordagem de penalização. 


21 OBJETIVO 


Em continuidade ao trabalho desenvolvido por (NASRI, 2010) e, também, em 
(FERREIRA, 2013), esta tese tem como objetivo principal propor um Método de La- 
grangeano Aumentado Exponencial para a resolução do Problema de Equilíbrio. Para 
tal, utilizaremos no lugar da penalidade clássica do método de lagrangeano aumentado 
proposto por (NASRI, 2010), a qual só possui derivada até primeira ordem, a penalidade 
exponencial, a qual possui derivadas de todas as ordens. 

Ainda, na tese de (FERREIRA, 2013) foram apresentadas duas abordagens 
de algoritmos do método de lagrangeano aumentado proposto em (NASRI, 2010). Tais 
algoritmos consistem na resolução de um sistema não linear e não diferenciável, originado 
do passo principal do algoritmo. A não a diferenciabilidade do sistema não linear é devido 
ao fato de que (NASRI, 2010) utlizam a penalidade quadrática clássica, conhecida na 
literatura, como PHR- Powell, Hestenes e Rockafellar. Por isso, a importância de se 
usar outras penalidades que possuam derivadas de ordem superior. É isso que estamos 
propondo, nesta pequisa, como uma das principais contribuições. 

Com a introdução da penalidade exponencial, o sistema não linear obtido pelo Passo 
2 do algoritmo de (NASRI, 2010), será diferenciável o que passa a ser um grande ganho 
em relação ao trabalho apresentado em (FERREIRA, 2013). Do ponto de vista teórico, 
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precisamos mostrar a equivalência entre lagrangeano aumentado exponencial e método 
de ponto proximal. 

Outras contribuições desta pesquisa são: mostrar que a sequência gerada pelo 
al- goritmo do novo método converge, implementar esta nova metodologia e testá-la em 
experimentos numéricos. 

Uma consequência inevitável quando se implementam algoritmos que utilizam 
a penalidade exponencial é o aparecimento de problemas de mau condicionamento. 
Para contornar esta possível dificuldade, utilizaremos uma estratégia de interpolação 
quadrática na implementação do algoritmo proposto. Por fim, para realizar experimentos 
numéricos, escolhemos na literatura duas classes de problemas: problemas de 
equilíbrio geral e problemas de equilíbrio de Nash generalizados, apresentados em 
(FACCHINEI; KANZOW, 2009). 


31 ORGANIZAÇÃO DO TRABALHO 


Neste primeiro capítulo apresentamos a motivação, objetivos e uma descrição da 
tese. 

No segundo capítulo apresentamos um apanhado geral sobre os métodos de 
penalidades, com destaque para sua forma geral e estrutura algorítmica. Inicialmente, 
destacamos os métodos de barreiras e, em seguida, os métodos de penalidades 
externas. O item principal desta parte é o método de lagrangeano aumentado, aqui 
apresentado para problemas de otimização. Ainda neste capítulo, dada a estreita 
relação entre lagrangeano aumentado e método de ponto proximal para otimização em 
Rê, descrevemos de forma resumida este importante método, com destaque para sua 
estrutura algorítmica e condições de convergência. O método de ponto proximal com 
distância de Bregman é destacado no final. 

Em seguida, no terceiro capítulo, considerando o trabalho de (NASRI, 2010), 
de onde retiramos o problema abordado nesta tese, nos voltamos para o problema de 
equilíbrio e, inicialmente, apresentamos o método de lagrangeano aumentado para 
este problema. Duas versões deste método são destacadas: extragaradiente inexato 
(IALEM) e ponto proximal extragradiente inexato (IPPEM). A versão ponto proximal 
extragradiente inexato irrestrito, como descrevem (NASRI, 2010), é apresentada com 
destaque para a fundamentação teórica que o acompanha. 

No quarto capítulo, descrevemos como principal contribuição deste trabalho: um mé- 
todo de lagrangeano aumentado exponencial aplicado ao problema de equilíbrio, seguindo, 
principalmente, o exposto em (NASRI, 2010) e (FERREIRA, 2013). Na composição desta 
nova metodologia, destacamos dois novos Algoritmos: IALEMe e IPPEMle. A prova de 
convergência é apresentada, assim como, a equivalência entre IALEMe e IPPEMle. A 
estratégia de penalização exponencial poderá acarretar problemas de mau condicionamento 
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devido ao rápido crescimento da função exponencial. Uma alternativa para contornar esse 
possível problema é fazer um ajuste quadrático da função penalidade. Esse procedimento 
é descrito na quarta seção deste capítulo. 

Com o intuito de validar a metodologia proposta, no quinto capítulo, realizamos 
experimentos numéricos e a análise dos resultados obtidos. Inicialmente aplicamos o Algo- 
ritmo IPPEMlIe para resolver problemas de equilíbrio de Nash Generalizado e problemas 
de equilíbrio geral, utilizando primeiramente a penalidade exponencial e num segundo 
momento utilizando a penalidade exponencial com ajuste quadrático. 


No último capítulo, apresentamos as considerações finais e conclusões. 
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LAGRANGEANO AUMENTADO PARA PROBLEMAS DE 
OTIMIZAÇÃO 


11 INTRODUÇÃO 


Os métodos de lagrangeano aumentado são baseados em penalização e utilizados, 
em geral, para resolver problemas de programação não linear com restrições. São métodos 
que buscam resolver o problema principal através de um processo iterativo de maneira 
que, a cada iteração, resolve-se um problema auxiliar irrestrito. Esse subproblema é 
obtido adicionando-se à função objetivo um múltiplo das restrições do problema original. 
Dentre as referências clássicas que tratam deste método, destacamos (POWELL, 1969; 
HESTENES, 1969; ROCKAFELLAR, 1974; BERTSEKAS, 1976; SANTOS, 1999; IUSEM, 
1995), entre outros. Uma descrição resumida de tais métodos será apresentada logo após 
a apresentação dos métodos de penalidades clássicos, os quais são descritos a seguir. 


21 MÉTODOS DE PENALIDADES 


Nesta seção, adotaremos como referências, considerando as diversas existentes, 
(SAN- TOS, 1998; LUENBERGER, 2008; GRIVA; NASH; SOFER, 2009), que apresentam 
conceitos, algoritmos e observações sobre os métodos de penalidades. A seguir 
apresentamos um resumo de tais métodos. 

Como tratado em (SANTOS, 1998), uma ideia bastante natural e frequentemente 
utilizada em programação não linear é tentar converter problemas complexos em outros, 
em geral mais simples, cuja resolução seja conhecida. Seria importante se um problema de 
otimização geral pudesse, de alguma forma, ser reduzido a um problema de minimização 
sem restrições, minimização em caixas ou politopos. 

Na chamada penalidade interna a função objetivo é modificada agregando um termo 
funcional que tende a infinito quando o ponto se aproxima da fronteira do conjunto viável. 
Forma-se assim uma espécie de barreira intransponível: métodos irrestritos começando 
no interior da região são desencorajados a se aproximar do contorno devido a valores 
muito altos da função objetivo. Por esse motivo, os métodos de penalidade interna são 
também conhecidos por métodos de barreira. Esses são, por outro lado, os mais antigos 
métodos de pontos interiores, com prestígio radicalmente incrementado após a revolução 
que seguiu-se ao trabalho de (KARMARKAR, 1984). 

Na penalidade externa, muitas vezes denominada simplesmente de penalidade, 
acrescenta-se à função objetivo um termo cujo custo aumenta com a violação das restrições. 
A solução de um problema penalizado externamente está, geralmente, fora do conjunto 
viável, mas se aproxima dele quando o termo de penalidade é muito grande. 

No trabalho de (SANTOS, 1999), os autores comentam que a razão pela qual a 
penalidade não é o procedimento universal para lidar com restrições é que o parâmetro 
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que deve multiplicar a função-restrição para castigar violação (na externa) ou o risco 
de violação (na interna) provoca, ao tomar valores extremos, mal-condicionamento do 
problema penalizado. Outra dificuldade que pode surgir ao adotarmos procedimentos de 
penalização é que a própria estrutura do problema é transformada quando uma restrição 
é acrescida à função objetivo, em geral, complicando a fisionomia desta. No entanto, 
todas as estratégias de penalidade estão vivas na otimização atual por sua simplicidade, 
adaptabilidade para problemas de grande porte. Esse detalhe é fortemente destacado em 
(FACCHINEI; KANZOW, 2009; FACCHINEI; KANZOW, 2010b), entre outros. 

Os métodos de penalidades utilizados para resolver problemas com restrições, 
são iterativos e cada iteração transformam o problema original em outro problema que 
é irrestrito. Eles servem tanto para problemas com restrição de igualdade como para 
restrições de desigualdade ou mesmo para aqueles que são formados por esses dois tipos 
de restrições. Por questões didáticas serão tratados dois tipos de métodos de penalidades: 
métodos de barreira e métodos de penalidade externa. 


2.1 Métodos de Barreiras 


Os métodos de barreira são métodos de penalidade que, em geral, são utilizados 
para resolver problemas com restrições de desigualdade da forma: 


Minimizar f(x) 


sa. 9X) <0, i=1,2,...,m (2.1) 
xe R” 
emquef:R"'>Reg:R"->Rparai=1,...,m, são as funções que formam o 
problema (2.1), sendo f a função objetivo e as funções 9, para i=1,...,m, as restrições 


que determinam o conjunto viável. 


A função penalizada em relação ao problema (2.1) é definida como 
(x, p)e R”ºx R++>—> B(x, 0) = f(x) + pp(x) (2.2) 
em que p é o parâmetro de penalidade e p : R” — R é uma função de penalidade, 
a qual é chamada, neste caso, de função barreira. Há dois tipos bastante conhecidos na 


literatura, a função barreira logarítmica e a função barreira inversa, as quais são definidas, 
respectivamente, por 


xeRv> p(x)=- > In (-g(x) (2.3) 


i=1 
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La | 
gi(m) 


xeRo> p(y)=-, 


i= 


(2.4) 


pa 


em que g, para i=1,...,m, são as restrições do problema (2.1). 
A seguir apresentamos o algoritmo de penalidade interna tipo barreira. 


Algoritmo 2.1 Algoritmo de Barreira Aplicado ao Problema (2.1) 
Dados pº >0, k=0 
Enquanto o critério de parada não for satisfeito, encontre: 

x“! e arg minfB(x, p9 | xe R7 

1 

Faça p**! < p* (por exemplo: p**! = 5) p*) 

k=k+1 
Continue. 


No Algoritmo (2.1), B é uma função barreira (2.2), que pode ser formada pela barreira 
logarítmica dada na relação (2.3) ou pela barreira inversa dada na relação (2.4). A tarefa 
mais pesada neste algoritmo é a resolução do subproblema gerado, o qual consiste na 
minimização da função barreira. No entanto, este subproblema é irrestrito e, desta forma, 
pode ser resolvido por qualquer método de programação não linear irrestrita, como por 
exemplo: Newton, Gradiente Conjugado, Região de Confiança, entre outros. 

Um ponto negativo dos métodos barreiras é que só se consegue provar a convergência 
se o parâmetro de penalidade p — 0. Isso poderá fazer com que os subproblemas gerados 
pelo método tornem-se mal condicionados, como vemos no exemplo a seguir. 

Considere o seguinte problema extraído de (SANTOS, 1998). 


Minimizar (x, + 1)2 + (x, — 1)? 


sa. -x,<0 
que tem solução x* = (0, 1)”. 


Neste caso considerando a função barreira p(x) = — In(x,) e a função penalizada 
(2.2), tem-se: 
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B(x, D)=(x,+1)2+(x,-1)2- pIn(x). (2.5) 
Para verificar o que foi dito anteriormente sobre o mau condicionamento da matriz 
Hessiana dos subproblemas gerados pelo Algoritmo 2.1, será calculado o gradiente e a 
Hessiana de B dado na relação (2.5). 


v.B(x, Pp) = Am +1)— p/m (2.6) 
P.B(x p)= e y (2.7) 


Em (2.6) os pontos estacionários, considerando x, >0 são da forma 
E ps 1) 
X'= , E 
2 


E, dado que p >0, tem-se em (2.7) que 


V2B(x', p) = p+i—v1+2p | : (2.8) 
0 2 


Como a convergência só é obtida fazendo p tender a zero, tem-se 


2p 
lim ——— | 
mo p+1- IFP 


E, desta forma, a matriz Vº?,.B(x*, p) em (2.8) é mal condicionada. Este foi um dos 


fatos que motivou o desenvolvimento de novos métodos de penalização chamados de 
lagrangeano aumentado. 


2.2 Métodos de Penalidades Externas 


Os métodos de penalizações externas são utilizados tanto para resolver problemas 
com restrições de desigualdade como para problemas com restrições de igualdade, ou 
seja, na forma 


min f(x) (2.9) 
sa. gx) <0ec(x)=0 


comg:Rr>R7mec:Rºr> Rº. 

A solução de um problema penalizado externamente está, geralmente, fora do 
conjunto viável, mas se aproxima dele quando o termo de penalidade é muito grande, 
(SANTOS, 1999). 

Se considerarmos problemas contendo somente restrições de igualdade, a 
penalidade utilizada com mais frequência na literatura é a quadrática, apresentada a seguir: 


Lagrangeano aumentado para problemas de otimização 


10 


1 
xeRo>p()=5 5 let? = 5 o(weo = 55 teor (2.10) 


e a função penalidade externa neste caso será da forma: 


p 
(x pe RixR++>— P(x, p)=fx)= 5 > (cg)? (2.11) 
t=1 


com p >0, fixado. Neste caso, p irá aumentar em todas as iterações. 
Por outro lado, para problemas com restrições de desigualdade uma penalidade 
utilizada é a seguinte: 


xeRo> p(x)= > max(0, g(x) (2.12) 
i=1 
Podemos observar que esta penalidade é não diferenciável em pontos do seu o 


domínio. Portanto, uma alternativa é substituir esta penalidade por: 


xeRo> p(x)= > (maxi0, 9, (0)? (2.13) 
1=1 


Neste caso, a função penalidade externa assume a forma: 


(gpeR'xR+>> P(x p)=f)+p >; (maxXO, 9, (x)? (2.14) 
Agora, considerando o problema (2.9)/Temos: 


m p 
(g peRxR+>—> P(x p)=HX)+p (tm gi)? + à Eleito?) (2.15) 
i=1 i= 
O algoritmo a seguir sintetiza o método de penalização externa aplicado ao problema 


(2.9). 

Algoritmo 2.2 Penalidade Externa Aplicada ao Problema (2.9) 
Dados: x e Rn, pº>0 

k=0 

Enquanto o critério de parada não for satisfeito, encontre: 

x" e arg min(P (x, p9) Ix e Rº 


Faça p'“! > p* (por exemplo: p**! = 20º) 
k=k+1 


Continue. 
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A penalidade P (x, p*) é dada em (2.15) e p* é o parâmetro de penalidade. 


31 LAGRANGEANO AUMENTADO 


Os métodos de lagrangeano aumentado foram inicialmente propostos, para 
problemas com restrições de igualdade, por (POWELL, 1969) e (HESTENES, 1969) em 
meados da década de 60. Foram chamados, nesta época, de métodos de multiplicadores 
e mais tarde, generalizados por (ROCKAFELLAR, 1974), para problemas com restrições 
de desigualdade. Descreveremos nesta seção os métodos de lagrangeano aumentado 
para problemas com restrições de igualdade e desigualdade. Tais métodos seguem a 
mesma filosofia dos métodos de penalidades. A diferença principal está na introdução dos 
multiplicadores de Lagrange, também conhecidos como variáveis duais. Desta forma, são 
geradas duas sequências, (x), chamada de primal e, £u), cnamada de dual. A sequência 
primal, em geral, não é viável (não satisfaz g(x) < 0), mas espera-se que vá se tornando 
viável à medida que k cresce. Já a sequência dual é sempre viável, no sentido da viabilidade 
dual, u* > 0, para todo k > 0. 

Segundo (MATIOLI, 2001), uma característica importante dos métodos de lagran- 
geano aumentado é a relação de dualidade que estes possuem com os métodos de ponto 
proximal. Afirma, ainda que, um dos grandes resultados alcançados durante a época em que 
estes métodos foram sendo resolvidos foi exatamente esta conexão entre as sequências 
primal e dual ou, em termos de métodos, do lagrangeano aumentado e de ponto proximal. 
Ainda segundo (MATIOLI, 2001), este resultado garante que, sob hipóteses adequadas, a 
sequência de multiplicadores gerada pelo algoritmo de lagrangeano aumentado é a mesma 
obtida por um certo algoritmo de ponto proximal aplicado ao problema dual. Destaca-se, 
ainda, que toda a teoria de convergência da sequência de multiplicadores gerada pelos 
algoritmos de lagrangeano aumentado é desenvolvida no dual, via métodos proximais e a 
parte de computação é feita toda no primal. 

Os métodos de ponto proximal fazem uma espécie de regularização somando à 
função objetivo um termo positivo, geralmente chamado de núcleo. Quando são aplicados 
ao problema dual do problema de programação convexa funcionam como uma espécie de 
barreira, pois forçam os pontos gerados a ficarem no ortante positivo. Veja uma descrição 
completa de tal método em (IUSEM, 1995). 

Nos métodos de lagrangeano aumentado mais modernos, desenvolvidos nos 
últimos anos, praticamente toda teoria de convergência é desenvolvida utilizando métodos 
de ponto proximal com quase-distâncias generalizadas tipo Bregman ou q-divergências. 
Por outro lado, a implementação destes métodos é feita no primal, isto é, os subproblemas 
resolvidos envolvem somente variáveis primais estando as variáveis duais fixas, como será 
visto adiante. 

Trataremos a seguir, em seções distintas, os métodos de lagrangeano aumentado 
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para problemas com restrições de igualdade e restrições de desigualdade. 


3.1 Lagrangeano Aumentado para Problemas com Restrições de Igualdade 


Nesta subseção apresentamos o método de lagrangeano aumentado em sua forma 
clássica, ou seja, em sua formulação semelhante àquela dada inicialmente por (POWELL, 
1969) e (HESTENES, 1969). 

O problema a ser tratado é da forma: 


Minimizar f(x) 


sa. clxy)=0, i=1,...,p (2.16) 
xe R” 
emque f:Rº'>Rec:R">Rparai=1,..., p, são funções continuamente 


diferenciáveis. 
A função lagrangeana associada ao problema (2.16) é definida por: 


(x A)e Rºx Rº>o> [(x, À) fx) + ATo() = f(x) + 3 AC,(x) (2.17) 
As condições necessárias de primeira ordem darihisia que se x* é uma solução 
ótima do problema (2.16), então existe um vetor de multiplicadores de Lagrange A* tal que 
V Mx, A*) = 0, em que I é a função lagrangeana (2.17). 
Desta forma, se fe c, para todo i=1,..., p, são diferenciáveis então os candidatos 
a pontos estacionários para o problema (2.16) devem satisfazer o seguinte sistema de n + 


p equações não lineares por n + p variáveis 


(2.18) 


Vf(x) +ATVe(x) = 0 
elsj=0 4=12. 0): 


Veremos que os multiplicadores de Lagrange no método de lagrangeano aumentado 
a ser apresentado nesta seção são atualizados de maneira que a primeira equação em 
(2.18) seja satisfeita. 


Definição 2.1 Dizemos que dois problemas A e B são equivalentes quando toda 
solução de A é também uma solução de B e vice-versa. 

No sentido da Definição 2.1 podemos verificar que o problema (2.16) é equivalente 
ao seguinte problema 


Minimizar f (x) + A” c(x) 
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sa. cXM)=0, i=1,...,p (2.19) 
xe R” 


Utilizando a penalidade quadrática apresentada em (2.10), tem-se a seguinte função 
penalizada para o problema (2.19) 
(x,M,p)€R"?xR xR, > L(x,),0) = f(x) + Ne(x) + Se" (x)e(x) = 


Sul (220) 


i—1 


to ]* 


= f(x) +D Meia) E 


da 


Neste caso, p é o parâmetro de penalidade e À é o vetor de multiplicadores de 
Lagrange associado às restrições de igualdade. Pode-se observar que na função (2.20) 
aparecem p variáveis a mais que na função penalidade externa (2.11) apresentada 
anteriormente neste capítulo. No entanto, com a abordagem proposta por (POWELL, 1969) 
e (HESTENES, 1969) estas variáveis são atualizadas por uma fórmula fechada. Assim, para 
determinar uma boa estimativa para os multiplicadores de Lagrange, deriva-se a função L 
dada em (2.20) com relação à variável x, estando À e p fixados, e iguala-se a zero, ou seja, 


VL(z,A,0) = V f(x) +) Ve(x) + pe” (2)Ve(x) = 
=Vf(x)+(N + pe(x)) Ve(x) = 0. (2.21) 


Desta forma, uma boa estimativa para os multiplicadores é 
A=A+ pe(x), (2.22) 


pois (2.21) teria a forma 


VE) +AVT (x) =0 (2.23) 


e esta última equação é exatamente a primeira equação dada na relação (2.18), 
ficando apenas a segunda equação, c(x) = O (c(x) = (c (x), ..., cl), para ser satisfeita. 

Logo, a metodologia proposta por Powell e Hestenes é a introdução de um algoritmo 
que gere candidatos à solução, no sentido de serem pontos estacionários do problema 
(2.16). Estes devem satisfazer o sistema (2.18). Portanto, atualizando-se os multiplica- 
dores de Lagrange pela fórmula (2.22) a primeira equação do sistema (2.18) é satisfeita 
automaticamente, enquanto a segunda é forçada a ser satisfeita através da atualização do 
parâmetro de penalidade. 

Posto isso, apresentamos o algoritmo de lagrangeano aumentado proposto por 
(POWELL, 1969) e (HESTENES, 1969). 
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Algoritmo 2.3 Lagrangeano Aumentado Aplicado ao Problema com Restri- ção 
de Igualdade 


Dados pº>0,hNºeR7,a>1,k=0 


Enquanto o critério de parada não for satisfeito, encontre: 


x" e arg min(L(x, AS, 09) | x e R'), com L dada em (2.20) (2.24) 


Faça: Nº! = Nº + p*e(x**") 


Se lex >0, 1le(x9Il 


Faça p*“*! = ap* (por exemplo: a = 2). 


Senão p“! = p* 


k=k+1 


Continue. 


Alguns comentários sobre o algoritmo: o subproblema (2.24), ou seja, dado M*e RP e 


p*>0, encontrar 


x“! e arg miníL(x, A, p9) : xe Rº) 


é chamado subproblema interno, que consiste na minimização de uma função 
irrestrita e a função L é o lagrangeano aumentado. 

A parte representativa de maior esforço computacional do algoritmo é encontrar 
uma solução para este subproblema, o que pode ser executado utilizando-se métodos de 
programação irrestrita, como Newton, Região de Confiança, dentre outros. 

Após a resolução do subproblema (2.24) atualizam-se os multiplicadores de 
Lagrange À, pela relação (2.22). Em seguida faz-se uma verificação da viabilidade, c(x) = 
0. Se houver ganho de viabilidade do ponto x* para o ponto x“! mantém-se o parâmetro de 
penalidade, caso contrário, o mesmo é aumentado. Continua-se o processo até que algum 
critério de parada seja satisfeito. 

Como o parâmetro de penalidade só é atualizado se as restrições estiverem “pouco 
satisfeitas”, o que não ocorre nos métodos de penalidades clássicos apresentados nas 
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seções precedentes, reduz-se a possibilidade de mau condicionamento dos subproblemas 
gerados pelo método, o que os torna vantajosos. 


3.2 Lagrangeano Aumentado para Problemas com Restrições de Desigualdade 


Na sequência será apresentado o método de lagrangeano aumentado proposto 
em (ROCKAFELLAR, 1973) (ver também (MINOUX, 1986) e (BERTSEKAS, 1976)), 
para o problema com restrições de desigualdade, generalizando desta forma o método 
de multiplicadores proposto por (POWELL, 1969) e (HESTENES, 1969), apresentado na 
subseção anterior. Além disso, serão apresentados métodos de lagrangeano aumentado 
mais modernos cuja teoria de convergência é feita utilizando-se o método de ponto proximal 
aplicado ao problema dual com alguma quase-distância (Bregman ou &-divergência), 
seguindo o exposto em, por exemplo, (IUSEM, 1995; IUSEM; TEBOULLE, 1995). 

O problema tratado nesta subseção é com restrições de desigualdade e tem a forma 


Minimizar f(x) 
sa. g9()<0, i=1,2,...,m, (2.25) 
xe R? 


em que, feg, parai=1,2,...,m, são funções definidas em R” com valores em R. 

Uma opção para a generalização do método de lagrangeano aumentado para o 
problema (2.25) é obtida através da introdução de variáveis de folga para as restrições de 
desigualdade g(x) = 0, onde (g(x) = (9,(X), 90%), ...., 949)”, conforme encontramos em 
(MINOUX, 1986; ROCKAFELLAR, 1973). Desta forma, o problema (2.25), com variáveis de 
folga s,>0, i=1,...,m, possui o formato: 


Minimizar f(x) 


sa. 9) +8s,=0, |=1,2,...,m. (2.26) 
s20,i=1,2,...,m 
xe R” 


Como foifeito na seção anterior, ao elaborar um algoritmo de lagrangeano aumentado 
aplicado ao problema (2.26), o subproblema gerado tem a forma: 


m 


Minimizar (rto — E jutadto) +s;) + 13 (gilz) +82 : 2xER",s> ) + A2BT) 


i1=1 2 1—1 
e, neste caso, yu é o vetor de multiplicadores de Lagrange e n é o parâmetro de 
penalidade. 
Para o parâmetro de penalidade positivo, mostra se que o problema (2.27) se reduz 
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à (ver (ROCKAFELLAR, 1974; MINOUX, 1986)): 


Minimizar (ta + rh P(gi(x), ui) : ZE | , (2.28) 
q—1 
e, neste caso, para i= 1,..., m, a função penalidade P é dada pela fórmula: 
(gi)? + mile) se gi(7) > E. 
Plo(r), om) =) —u2 y (2.29) 


Cc. 


4m 


Note que a variável de folga s não aparece em (2.28) e nem em (2.29). Portanto, 
o algoritmo de lagrangeano aumentado desenvolvido para o problema com restrições de 
desigualdade é análogo ao Algoritmo 2.3, porém com subproblema interno dado por (2.28) 
e (2.29). Desta forma, substituindo a penalidade P dada em (2.29) na função (2.28) temos 
a função lagrangeano aumentada proposta por Rockafellar 


m 


1 
(x,u,n) ER" x R$ xRy > L(z,u,n) = f(x) + no [(max(0, ui + 2ngi(z))? — pê] 
= 
(2.30) 


Algoritmo 2.4 Lagrangeano Aumentado para o Problema com Restrição de 
Desigualdade 
Dados yº e R”, (n,) uma sequência de termos positivos e k = 0. 


Enquanto o critério de parada não for satisfeito, encontre: 


x“ = arg min(L(x, uº, n,) | xe R') com Ldada em (2.30)) (2.31) 
Faça: 
ut = max(0, ué + 2n«g(X+1), parai=1,...,m 
k=k+1 
Continue. 


3.3 O MÉTODO DE PONTO PROXIMAL PARA OTIMIZAÇÃO EM R” 

Na seção anterior descrevemos um método de lagrangeano aumentado aplicado ao 
problema de otimização com restrições de desigualdade, utilizando variáveis de folga para 
transformar as restrições do problema (2.25) em igualdades. Na sequência, descrevemos 
outra metodologia baseada nos métodos de ponto proximal, conforme descrito em (IUSEM, 
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1995). 

Aideiaé baseada nos métodos primais-duais, ou seja, apresentaremos a convergência 
do método de ponto proximal para um problema irrestrito e depois, a equivalência entre as 
sequências duais geradas pelos métodos de ponto proximal e lagrangeano aumentado. 

Primeiramente, descrevemos o método de ponto proximal aplicado ao seguinte 
problema irrestrito 


Minimizar f(x) 
sa:xeR” (2.32) 


no qual f: R” — R é uma função convexa e limitada por baixo. 

Mesmo o problema sendo convexo e f limitada por baixo, não existe a garantia de 
um único minimizador. Para assegurar a unicidade do minimizador é acrescentado à função 
objetivo do problema (2.32) um termo de regularização. Assim, o método de ponto proximal 
gera uma sequência (x da seguinte maneira: 


ei (2.33) 
x*'= argmin, af (09) + a llx— x 2) (2.34) 


em que a, é um número real satisfazendo O <a, < a para algum a >0, incluindo o 
caso a, constante. 

O resultado a seguir estabelece a convergência da sequência (x, gerada pelo 
método de ponto proximal (2.33) e (2.34). 


Teorema 2.1 Convergência do método de ponto proximal 

Considere f: R" > R uma função convexa e continuamente diferenciável. Assuma 
que o conjunto solução U de minimizadores de fem R" é não vazio. Então a sequência 
gerada por (2.33) e (2.34) converge para um ponto xe U. 


Prova. Ver (IUSEM, 1995). 

Para provar a convergência da sequência dual (u'“, gerada pelo Algoritmo 2.4, será 
necessário introduzir o problema dual associado ao problema primal (2.25). Para isso, 
considere fe g,convexas e diferenciáveis, para todo i= 1,2..., m. Associado ao problema 
(2.25) define-se a função lagrangeana clássica /: Rr x Rr — R por 


f(x) +> malz) seu>0 (2.35) 
i=1 


+00 G:6; 
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Considere a função objetivo W : Rm e R definida por 


yu) = min Mx, py) (2.36) 
xeR” 


em que /é a função lagrangeana (2.35). 

O resultado seguinte prova que a sequência (u'), gerada pelo Algoritmo 2.4, é a 
mesma sequência gerada pelo algoritmo de ponto proximal (2.33) e (2.34) quando aplicado 
ao problema —wy. 


Teorema 2.2 Equivalência entre ponto proximal e lagrangeano aumentado - para 
problema de otimização 

Considere (uy a sequência gerada pelo algoritmo de ponto proximal, (2.33) e 
(2.34), aplicado ao problema Aim — Wu), coma, = as tu” a sequência dual gerada pelo 


nk 
algoritmo de lagrangeano aumentado 2.4. Seyº = yº então “= y* para todo k. 


Prova. Ver (IUSEM, 1995). 

Até este momento poderíamos nos dar por satisfeitos, pois como mencionado nos 
parágrafos precedentes, os métodos de lagrangeano aumentado para os problemas com 
restrições de igualdade e desigualdade apresentados até aqui, contornam, de certa forma, 
o problema do mau condicionamento gerado pelos métodos de penalização clássicos. 

Com a introdução dos multiplicadores, ou variáveis duais, foi possível mostrar a 
convergência (dual no caso de desigualdade) dos métodos de lagrangeano aumentado, 
sem que o parâmetro de penalidade aumentasse em todas as iterações. Inclusive, este 
parâmetro pode ser constante para o método de lagrangeano aumentado, proposto por 
(ROCKAFELLAR, 1974), para problemas com restrição de desigualdade como pode ser 
visto no Algoritmo 2.4. 

A questão principal está na determinação da solução do subproblema gerado pelo 
Algoritmo 2.4, ou seja, na solução do subproblema dado na relação (2.30). Note que 
esta função só possui derivada, em todo o seu domínio de definição, até primeira ordem. 
Portanto, métodos que usam informação de segunda ordem, como é o caso do método de 
Newton e de região de confiança, não podem ser utilizados, ou se forem utilizados podem 
gerar instabilidade quando os pontos obtidos se aproximarem de valores em que a função 
não é diferenciável. Desta forma, a comunidade de otimização desenvolveu métodos em 
que as funções envolvidas possuam derivadas de ordem superiores. 

O desenvolvimento seguinte visa explorar a relação primal-dual. Mostraremos a 
convergência do método de ponto proximal e, em seguida, a equivalência com o método de 
lagrangeano aumentado. Porém agora, ao invés da regularização quadrática apresentada 
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na função (2.34) do método de ponto proximal, utilizaremos uma quase distância mais 
geral. Da mesma forma, no primal deveremos substituir a penalidade quadrática da função 
lagrangeano aumentado (2.30) por uma quase distância que possua derivada de ordens 
superiores. 

As quase-distâncias mais utilizadas são: distâncias de Bregman e q&-divergência. 
Nesta tese, em função de sua especifidade, será apresentado somente o caso de distâncias 
de Bregman. Para isso, o desenvolvimento apresentado está baseado em (IUSEM, 1995). 

Diferentemente da regularização utilizada no método de ponto proximal quadrático, 
aqui será necessário que este termo penalize pontos inviáveis, pois o interesse é aplicar o 
método de ponto proximal ao problema dual do problema (2.25), no qual as variáveis duais 
são não negativas. Sendo assim, a regularização funcionará como nos métodos de barreira, 
descritos no início deste capítulo. O objetivo da regularização é penalizar “infinitamente” 
pontos que se aproximem da fronteira. Para isso, considere o seguinte problema restrito: 


Minimizar f(x) 
sa xeS (2.37) 


em que Sc R"é fechado. 


Seguindo o exposto em (IUSEM, 1995), uma abordagem para resolver (2.37) é a 
introdução de uma função de penalidade p satisfazendo p(x) = sex é S. Paraa >0,0 
problema 


min Lf (09) + ap(x)) (2.38) 
xeR” 


terá sua solução em Ss (se as soluções existirem), como foi o caso desenvolvido 
para o método de ponto proximal clássico. Se (2.37) tiver solução e p for apropriadamente 
escolhida, então (2.38) terá uma única solução x(a) e sob hipóteses adicionais pode ser 
provado que: lim x(a) existe e resolve (2.37). 


a-o 

Como no caso da regularização (aqui olhando também como penalização clássica), 
o problema é que (2.38) se torna mal condicionado quando a se torna muito pequeno, 
isto é, para pontos próximos da fronteira de S, p(x) se torna muito grande e ap(x) se torna 
indeterminado quando a se aproxima de 0. Logo, a ideia é fazer uma combinação de 
penalização com a abordagem de ponto proximal para obter convergência para a distante 
de zero. Para alcançar este objetivo, usaremos uma função quase-distância D: Sx S — 


R+ (em que Sé o interior de S) tal que D(x, y) = 0 se x= ye se aproxima de infinito quando 
y se aproxima da fronteira OS de S. 
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3.4 Funções e Distâncias de Bregman 


Considere S c R" não vazio, aberto, convexo e seu fecho S, d uma função convexa 
realemSe D,:Sx S-— R definida por 


Dix, y) = dx) — dy) — Va)” (x — 3). (2.39) 


A função d será chamada uma função de Bregman e D uma distância de Bregman 
induzida por d, se cumprir as seguintes condições: 

B1) dé continuamente diferenciável em S. 

B2) d é estritamente convexa e contínua em S. 

B3) Para ô e R, os conjuntos de nível parciais 


(x 0)=(xe S: Dix )<ôerdxô)=(yeS:D/x y<ô), 


são limitados para todo y e Se todo x e S, respectivamente. 

B4) Se (y9% c S convegir para y* então D (y*, y*) converge para 0. 

B5) Se (x) c Se (y) c S são sequências tais que (x é limitada, lim y*= y* e 

lim Da0&, y)= 0 então lim = y. pi 

De acordo com (IUSEM, 1995), vale a pena destacar algumas observações a 
respeito da definição e propriedades anteriores: 

(i) O conjunto S é chamado de zona de d. 

(ii) D(x, y) O para todo x e Seye S. Além disso, D(x, y) = 0 se, e somente se, x 
— y. 

(ii) As propriedades B4 e B5 precisam ser verificadas somente nos pontos da 
fronteira OS de S, pois quando x*e y* estão em S elas se verificam automaticamente, como 
uma consequência de B1, B2 e BS. 

Uma função de Bregman d será chamada de coerciva na fronteira se a seguinte 
propriedade for satisfeita: 


B6) Se(y94 c Sétal que lim y*= y cos, então lim Vd(y9” (x — y*) = —oo para 

todo xe S. lã há 

A seguir são apresentados alguns exemplos de funções de Bregman que podem ser 
encontrados em (IUSEM, 1995). 

Exemplo 2.1 S = Rº, d(x) = x'Mx com Me R”” simétrica e positiva definida. Neste 


caso, Dix, )=(x-y'M(x-y)=Ilx- yll2, (Il, é a norma induzida por M). 
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n 


d(x) = Anos x, extendida com continuidade em oR 1 A 
= mn 
e convenção que O log O = 0. Neste” caso, DÁx, y) = Ds Jog É +yj — Z;), a qual é 
g=. É] 
conhecida na literatura como entropia (ou divergência) de Kullback-Leibler. 


Exemplo 2.22 S=R” 


++ 


n 


Exemplo 2.3 S — Ri,, d(x) = >o(zg — 28) coma >1,0<8B<1. Paraa-2e 
j=1 


ERR vm) eparaa=1eB=1/2 


Da(x,y) = 5 


= Ev vis. 


3.5 O Método de Ponto Proximal com Distância de Bregman 


Considere S c R” aberto e convexo, So fecho de Se f: R" — R umafunção convexa 
e contínua em S. O método de ponto proximal aplicado ao problema (2.37) é definido por: 


Xes (2.40) 
x"! = argmin, df (x) + a,D (x, x9) (2.41) 


em que d é uma função de Bregman com zona Se 0 <a, <a para algum a >0. 
A seguir trataremos da convergência do método de ponto proximal. Para isso assu- 
mimos que a função f que define o problema (2.37) é limitada por baixo em S. 


Teorema 2.3 Convergência do método de ponto proximal com distância de Bregman. 
Se o problema (2.37) tem soluções e d é coerciva na fronteira com respeito a S, então a 
sequência (x*) gerada por (2.40) e (2.41) converge para uma solução x* do problema (2.37). 

Prova. Ver (IUSEM, 1995). 

Apesar da existência de uma gama de metodologias de lagraneano aumentado 
utilizando distâncias de Bregman, aqui trataremos somente o caso do método exponencial 
apresentado por Bertsekas conforme (IUSEM, 1995) (ver também (BERTSEKAS, 1976)), 
o qual será de grande valia para nos auxiliar no desenvolvimento da nossa proposta, a ser 
apresentada no próximo capítulo. 

Considere novamente o problema de otimização com restrições de desigualdade, 
ou seja 


Minimizar f(x) 
sa. 9X) <0, i=1,2,...,m, (2.42) 


xe Rº 
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no qual fe g, parai=1,2...,m, são funções convexas e diferenciáveis. 
A função lagrangeano aumentado com penalidade exponencial, associada ao 
problema (2.42), é dada por: 


m 


(gumeRxRExR o Lou m=t0)+> Eu eNsite) (2.43) 


onde n > 0 é o parâmetro de penalidade e y é o vetor dos multiplicadores de 
Lagrange. 

Note que a função lagrangeano aumentado exponencial (2.43) é tantas vezes dife- 
renciável quanto forem fe g. O mesmo não se verifica no caso da função lagrangeano 
aumentado quadrática de Rockafellar considerada na relação (2.30). 

Fixada uma sequência de termos não negativos (n,), o método de lagrangeano 
aumentado exponencial gera as sequências (x e (uy, através dos seguintes passos: 


vºeRTeparak=0,1,2,... (2.44) 
x" = argmin, a 4x, vu Sn) (2.45) 
= uee(x), [=1,2,...,M (2.46) 


em que L é a função lagrangeano aumentado exponencial (2.43). 

Equivalentemente ao Teorema 2.2, é possível mostrar a convergência da sequência 
dual, 4u*%, gerada pelo método de lagrangeano aumentado exponencial, através da 
equivalência entre ponto proximal com distância de Bregman e lagrangeano aumentado. 


Para isso, recorremos ao seguinte resultado: 


Teorema 2.4 Teorema de equivalência para o problema de otimização com Bregman 
Considere 4u'*) a sequência gerada por (2.44)-(2.46) e (uk) a sequência obtida através de 


a 1 
nº = argmin,zot=U(p) + — Da(pt, pn) 
Mk (2.47) 


com y(u) = mintf(a )+ Ema )jed:R” > R dada por d(x = Suyloga;. Seu? 
j=1 
=yº então u* = y* para todo k. 


Prova. Ver (IUSEM, 1995). 
Note que, na equação (2.47), foi necessário escolher o parâmetro de regularização 
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ak = + sendo n, o parâmetro de penalidade no método de lagrangeano aumentado 
n 


exponencial. Esta escolha é necessária para se conseguir demonstrar tal equivalência. 
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LAGRANGEANO AUMENTADO PARA PROBLEMAS DE 
EQUILÍBRIO 


11 INTRODUÇÃO 


Atualmente o problema de equilíbrio tem recebido muita atenção dos pesquisado- 
res, comprovado pelos vários resultados teóricos e experimentos publicados. Entre estes 
estudos destacamos: (BIANCHI; PINI, 2005; BIGI et al., 2013; BLUM; OETTLI, 1994; 
BURACHIK, 2012; FACCHINEI, 2003; FLâM; ANTIPIN, 1996; IUSEM; NASRI, 2007; 
NASRI, 2010; IUSEM; KASSAY; SOSA, 2009; IUSEM; SOSA, 2003a; IUSEM; SOSA, 
2003b; MASTROENI, 2000; MATIOLI; SOSA; YUAN, 2012; MUU; OETTLI, 1992; MUU; 
QUOC, 2009; FACCHINEI; KANZOW, 2007; QUOC; MUU, 2004; TRAN; DUNG; NGUYEN, 
2008; SANTOS; SCHEIMBERG, 2011; NGUYEN; STRODIOT; NGUYEN, 2009; NGUYEN; 
STRODIOT et al., 2009). Em particular (BIGI et al., 2013), faz uma descrição da situação 
atual no que se refere a algoritmos e condições de existência de soluções. 

O problema de equilíbrio fornece uma abordagem unificadora de vários problemas 
matemáticos. Entre estes, inclui como casos particulares, problemas de otimização convexa, 
problemas de ponto fixo, problemas de complementaridade, problemas de equilíbrio de 
Nash, problemas de desigualdade variacional e problemas de minimização vetorial, como 
destacado em (BIGI et al., 2013; BLUM; OETTLI, 1994; IUSEM; SOSA, 2003b). 

Neste capítulo apresentaremos o problema de equilíbrio e uma abordagem recente 
apresentada por (NASRI, 2010), sobre o método de lagrangeano aumentado para a sua 
resolução. Neste artigo, os autores apresentam métodos de lagrangeano aumentado para 
resolução de problemas de equilíbrio em dimensão finita, cujos conjuntos viáveis são 
definidos por inequações convexas, generalizando o método de lagrangeano aumentado 
para otimização restrita. A cada iteração as variáveis primais são atualizadas resolvendo-se 
um problema de equilíbrio irrestrito, enquanto as variáveis duais são atualizadas por uma 
fórmula fechada. Além disso, é realizada a análise de convergência do algoritmo. 

O problema de equilíbrio com as propriedades mínimas necessárias para o desenvol- 
vimento teórico é apresentado na definição a seguir. 


Definição 3.1 Considere K um subconjunto de R” não vazio, convexo e fechado e 
uma função f: Kx K— R satisfazendo 

Pi: f(x,x)=0 para todo xe K, 

P2: f(x, ): K> Ré convexa e semicontínua inferior para todo xe K, 

P3:f(,): K> Ré semicontínua superior para todo ye K. 


O problema de equilíbrio, denotado por EP (f, K), consiste em determinar x e Ktal 
que f(x*,y) >0 para todo ye K. 
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O conjunto de soluções para EP (f, K) será denotado por S(f, K). 

Em termos de métodos computacionais para problemas de equilíbrio, (BIGI et 
al., 2013; FACCHINEI; KANZOW, 2010a) comentam que existem poucas referências 
encontradas na literatura. Este foi um dos fatos que serviu de motivação para buscarmos o 
desenvolvimento de algoritmos para resolução destes problemas. 

Neste cenário, destacamos o artigo de (NASRI, 2010), no qual propõem versões 
exatas e inexatas de métodos de lagrangeano aumentado para resolução do EP (f, K) em 
R", para o caso nos quais o conjunto viável K é da forma: 


K=(xe Rº:h(9) <0 (Isism) (3.1) 


em que h,: R” —» R são funções convexas para todo = 1,...,m. 

Observe, na definição do problema de equilíbrio, a dependência das duas variáveis, 
xe y, estarem no conjunto convexo K, ou seja, quer se encontrar x* que esteja em Ktal que 
f(x*, y) > 0, para todo y em K. Devido ao fato das variáveis estarem no conjunto convexo 
K, a maioria dos artigos encontrados na literatura acabam fazendo algum tipo de projeção. 
É bem conhecido que fazer projeção reduz-se em outro problema de otimização que pode 
ser de difícil solução, a menos que o conjunto K seja tipo poliedros (descritos por restrições 
lineares ou limites das variáveis). Com a introdução do método de lagrangeano aumentado 
esta dependência pode ser enfraquecida através da abordagem de penalização. 
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Assumiremos, como em (NASRI, 2010), que o conjunto K que define o problema de 
equilíbrio, satisfaz a Condição de Qualificação de Slater dada na definição a seguir. 


Definição 3.2 Condição de Qualificação de Slater 

Dizemos que o conjunto K, dado em (3.1), satisfaz a Condição de Qualificação de 
Slater se existir x e R”, tal que h(x) O paraie le h(x) <O para i É |, sendo Io conjunto de 
índices (possivelmente vazio), tal que a função h, é afim. 


Vale ressaltar que existem outras condições de qualificação que poderiam ser 
utilizadas. 

Para maiores detalhes consultar, por exemplo, (EUSTÁQUIO, 2007). 

A bifunção lagrangeana L : (R” x R7”) x (R" x R7) > R, associada ao problema de 
equilíbrio EP (f; K), (3.1), é definida por: 


Ele) (8,80) = 12,9) +35 Mi) — DO a (32) 


i=1 
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Destacamos que o problema de otimização (2.42) pode ser tratado como um caso 
particular do problema de equilíbrio EP (f, K), tomando-se f (x, y) = f(y) — f(x). Nesse 
caso, (3.2) reduz-se à L((x, À), (Y, //)) = Ky, À) — Kx, 4), sendo / a função lagrangeana (2.35) 
associada ao problema (2.42). 

Na sequência, apresentaremos o método de lagrangeano aumentado proximal 
para EP (f, K), conforme apresentado em (NASRI, 2010). Considere K o conjunto dado na 
definição (3.1), a função penalidade p,: R"x R"x R7 x R, — R definida por: 


2 2 
Pi(z,y,A,Y) = > (ms (o A+ ao) — (nas o A+ o) | (3.3) 


e a função lagrangeano aumentado proximal L:R"x R"xR”x RºxR, — R definida 
por: 


L(a, y, A, z, 7) = fiz, y) E yta Seg) 2) Tr Dpi(z, y, A, 7). (3.4) 


i=1 


Note que a função lagrangeano aumentado (3.4) torna o problema de equilíbrio EP 
(£ K) irrestrito, ou seja, as restrições que descrevem o conjunto K são incorporadas nesta. 
Desta forma, determina-se uma solução para o problema de equilíbrio restrito através 
da resolução de uma sequência de problemas irrestritos, utilizando-se a metodologia 
delagrangeano aumentado proximal. 

A seguir descrevemos os dois principais algoritmos apresentados em (NASRI, 2010) 
para o problema de equilíbrio EP (f, K): lagrangeano aumentado ex- tragradiente inexato e 
ponto proximal extragradiente inexato. 


Algoritmo 3.1 Lagrangeano Aumentado Extragradiente Inexato (IALEM) 


Passo 1: Considere uma sequência limitada ty ) c R... (X,A) e R" x R”, uma 


++? 


tolerânciad e (0,1) ek=0. 


Passo 2: Encontrar um par (X*, e) e R" x Rº tal que X* resolve EP (LSR?) com (LS, 
definida como: 


Lilz,y) = f(x,y) a vela o 2, y 2) + Sopi(z, y, a Ye) E (ey nai $), 


i=1 
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com p, dada por (3.3) e e* dado por: 
el < Oy 1X! — xa. 


Passo 3: Atualize a variável dual 


ha(T 
AH = max (o o eg E die) si Ldam): 
Y 
Passo 4: Se (x*, Nº = (X*, At!) pare. Senão faça: 
l 
xt+1 = fes —es, 
Vk 


k=k+1 evolte ao Passo 2. 


Note que no Passo 2 do Algoritmo IALEM, a função lagrangeano aumentado regula- 
rizada apresenta o termo (e”, y — x), O qual representa o erro. Daí a presença do termo 
inexato na denominação do algoritmo. 

A seguir apresentamos o Algoritmo de Ponto Proximal Extragradiente Inexato, o qual 


será utilizado como resultado auxiliar. 
Algoritmo 3.2 Ponto Proximal Extragradiente Inexato (IPPEM) 


Passo 1: Considere uma sequência limitada ty ) c R,., Xe K, uma tolerância 
de(0,1)ek=o0. 


Passo 2: Encontre um par (X*, e) e R" x Rº tal que X* resolve EP (f é K) com 
Ex p=t(xgp+yax- xs y-M)-(esy-x elledl < Oy kt xa. 


Passo 3: Se X* = x*, então pare. Senão faça: 
l 

xt = sé ms =sê", 
Vk 


k=k+1, 
e volte ao Passo 2. 


A análise de convergência do Algoritmo IPPEM é realizada usando as técnicas já 
conhecidas dos métodos de ponto proximal. Desta forma, (NASRI, 2010) mostram que o 
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Algoritmo IPPEM converge para um ponto x* e S(f, K). 

Já a prova de convergência do Algoritmo IALEM é realizada também pelos procedi- 
mentos padrões conhecidos na literatura de programação não linear, através da equivalência 
entre primal-dual. Para isso, são necessárias duas hipóteses adicionais: 


P4: Considere a função f dada em EP (f, K), 8 sub monótona, isto é, existe O > 0 tal 
que f(x,y) + f(y ) <0lx- yl2, vx, ye Rº. 


q 
P5: Paratodo xX),..., xe Rretodot,..., EA. tal que tel, se cumpre que 


1=1 
ui (= = tuo) < 0. 
i=1 k=1 
No próximo teorema apresentamos a prova de convergência do Algoritmo IPPEM. 
Teorema 3.1 Teorema de convergência do IPPEM 
Considere EP (f, K) satisfazendo P1-P5 e uma sequência exógena ly ) c (0, Y] para 
algum Y > 6, sendo 6 a constante de monotonicidade em P4. Se EP (f, K) tem solução e 


(x) é a sequência gerada pelo Algoritmo IPPEM, então (x converge para alguma solução 
x* de EP (£ K). 


Prova. Ver Teorema 5.8 de (IUSEM; NASRI, 2007). 
No Teorema 3.1, o problema de equilíbrio considerado é EP (f, K) significando a 
dependência do conjunto K. Na próxima proposição esta dependência será relaxada, ou 


seja, o Algoritmo IPPEM será aplicado para resolver EP (L, R" x R7), com L dada em (3.2). 


Proposição 3.1 Considere f satistazendo P1-P5 em R" xR" e Kdado por(3.1). Então 
L, como definida em (3.2), satisfaz P1-P5 em (R" x R”) x (R” x R7). 


Prova. Ver (NASRI, 2010) 
No Passo 2 do Algoritmo IPPEM será necessário resolver o seguinte subproblema: 


Obter um par (x*, e) e R" x R"tal que x K resolve EP (f,, K) sendo 
EO n)=t0))+yax- xs y-x) — (es y— >), (3.5) 
e o erro satisfazendo 


els dy llx* — xl. (3.6) 
Poderemos, agora, aplicar o Algoritmo IPPEM, para resolver o problema de equilí- 
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brio irrestrito EP (L, Rº x R7). Para isso, na iteração k deste algoritmo será necessário 
determinar a função de regularização L deste novo problema. Utilizando a relação (3.5) 
com L no lugar de ftemos: 


Lot A), (y 1)) = 

= L(x, A), (4, 4)) + val À) — (x AY, (4 40) — (x, A)) — (e% (y, 40) — (x, A) = 

=L((x, À), (4 2)) + yax— x AA + y4y— x >A) — (es y— x). (3.7) 

em que L é dada por (3.2) e não foi usado o erro associado a À e // no argumento de 
L. Isto está relacionado ao fato que no Passo 3 do Algoritmo IALEM a variável À é atualizada 
por uma fórmula fechada e, desta forma, é assumido que a atualização é realizada de uma 
maneira exata. 

O Algoritmo IPPEM para resolver EP (L, R” x R”) será denotado daqui em diante 
por IPPEMI, ou seja, Ponto Proximal Extragradiente Inexato Irrestrito, o qual é dado por: 


Algoritmo 3.3 Ponto Proximal Extragradiente Inexato Irrestrito (IPPEMI) 


Passo 1: Considere uma sequência limitada ty j c R,., (e, À9) eRrxR7, e R, 
k=0 e uma tolerância à e (0, 1). 


Passo 2: Encontre (8, Às eRrxRme(e40)eR"xR” tal que (8, NO) resolve 
EP(Le, Rºx R”) com 


Le((x, A), (4,20) = LM, À), (4,40) + VAX XY) + VMA AS UA) (E y 3). 


els dy (E — x, = AY 


A A ' = 
Passo 3: Se x*= x*e A*= A*, então pare. Caso contrário faça 


1 E 
qt = ed e nt (3.8) 
O próximo resultado, apesar de tratar da convergência do Algoritmo IPPEMI, foi 
apresentado em (NASRI, 2010) como um corolário, em razão do Algoritmo IPPEMI ser uma 


decorrência do Algoritmo IPPEM. Lembrando que a função Le e do Algoritmo IPPEMI, dada 
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na relação (3.7), foi obtida da função É do Algoritmo IPPEM, dada na relação (3.5). 


Corolário 3.1 Convergência do Algoritmo IPPEMI 

Considere EP (f, K) com K dado por (3.1), f satisfazendo P1-P5 em Rº x R" e uma 
sequência ty,» c (0, Y ] para algum Yy >6, em que 8 é a constante de sub monoionicidade 
de f. Seja £(x*, A9) a sequência gerada pelo Algoritmo IPPEMI. Se o problema EP (L, R" x 
R”) tem solução, enião ((x*, A) converge para algum par (x*, Nº) e S(L, Rº x R7). 


Prova. Ver (NASRI, 2010). 

Antes de tratar da equivalência da sequência dual gerada pelos Algoritmos IPPEMI 
e IALEM será necessário estabelecer a associação entre o conjunto de soluções de EP (f, 
K)e de EP(L, R" x RT). Esta relação será útil para garantir a unicidade de solução no passo 
principal de IPPEMI. Para isso serão necessários alguns resultados preliminares, os quais 
descrevemos na sequência. 

Para cada x e R”, define-se F.: R” — R como 


EO) = f(x 9). (3.9) 


a qual é convexa em decorrência da propriedade P2. Esta função será usada para 
estabelecer a relação entre S(f, K) e S(L, R" x R7). 


Proposição 3.2 Considere EP (f, K). As duas afirmações seguintes são equivalentes: 


Dx e SL KH. 
ii) x* minimiza F,. em K, sendo F,. dada por (3.9). 


Prova. Suponha que x* e S(f, K). De (3.9) e P1 tem se 
Fo) = 06, y) == [(00,%) = F, (1) 


JE 


para todo y e K, estabelecendo (i). Agora suponha que (ii) é satifeita. Usando 


novamente P1 e (3.9), tem se 
FO, p=) =f(X,x)=0 
para todo y e K, o que fornece o resultado desejado. 


Para os próximos resultados será necessário o conceito de par ótimo para o problema 
EP (f, K), ou seja: 
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Definição 3.3 O par (x*, Nº) e R" x R” será chamado um par ótimo para EP (f, K) se 


0 € 9Fre(1*) + DE, MOhi(a*), (3.10) 
N'>0 (Isism), (3.11) 
h(X)<0 (I<ism), (3.12) 
Nh(X)=0 (1 <ism), (3.13) 


em que OF, (x'), dh(x*) denotam o subdiferencial das funções convexas F. e 
h, respectivamente, no ponto x*. Observe que (3.10) - (3.13) são as condições de KKT 
associadas ao problema 


Minimizar Fx* (x) sujeito ax e K, 


sendo K dado por (3.1). 
As próximas proposições são resultados auxiliares para se obter a equivalência 
entre as soluções de EP (f, K) e de EP(L, Rº x RM). 


Proposição 3.3 Considere EP (f, K) e assuma que fcumpre P1-P3 em R" xR”. Então 
as seguintes afirmações são equivalentes. 

1) (xX*, A*) é um par ótimo de EP (f, K). 

ii) (6, A) e SL, Rº x R7). 

Prova. Ver (NASRI, 2010) 


Proposição 3.4 Considere EP (f, K) e assuma que f cumpre Pi-P3 em Rº x Rº. Se 
x* minimiza F.. em K, com F.. como em (3.9) e a condição de Qualificação de Slater dada 
na definição 3.2 se verifica para as funções h, que definem o conjunto K, então existe N* e 
R” tal que (x*, A*) é um parótimo para EP (f, K). Reciprocamente, se (x*, A*) é um par ótimo 
para EP (f, K), então x*minimiza F,. em K, com F.. como em (3.9). 

Prova. Ver (NASRI, 2010) 


Corolário 3.2 Considere EP (f, K) e assuma que f cumpre P1-P3 em R" x Rº. Se 
(X,A) e S(L, Rº x R”), então x* e S(f, K). Reciprocamente, se x* e S(f, K) e a Condição de 
Qualificação de Slater dada na definição 3.2 se verifica, então existe N* e R” tal que (x*, 
Ne S(L, Rº x Rº). 


Prova. Ver (NASRI, 2010) 
O Corolário 3.2 mostra que resolver EP (L, R" x R”) é suficiente para resolver EP 
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(£ K). Ou seja, este resultado garante que resolver o problema de equilíbrio irrestrito é 
suficiente para se resolver o problema de equilíbrio restrito. 

O resultado mais importante deste capítulo será o teorema de equivalência, ou seja, 
o teorema que garante que a sequência gerada pelo Algoritmo IPPEMI para resolver EP 
(L, R” x R7) coincide com a sequência gerada pelo Algoritmo IALEM para resolver EP (f, 
K). Antes, porém, apresentaremos um resultado técnico que será utilizado na prova de 
equivalência. 


Proposição 3.5 Considere f satisfazendo P1-P4, ze Rrey >0, sendo 6 a constante 
de sub monotonicidade de f dada em P4. Sef : Kx K> R é definida como 


Fx p=tixp+yx-zy-m-(ey-», 
então EP (F. , K) tem uma única solução. 


Prova. Ver (NASRI, 2010) 

Com o que foi desenvolvido até aqui tem se todas as ferramentas necessárias 
para apresentar o teorema de equivalência. Apesar deste teorema estar demonstrado em 
(NASRI, 2010), exibiremos sua prova aqui também, por se tratar de um resultado importante 
na abordagem que estamos apresentado nesta tese. 


Teorema 3.2 Equivalência para o Problema de Equilíbrio 

Considere EP (f, K) satisfazendo as propriedades P1-P4 em R" x R", uma sequência 

(yo c R,, uma tolerância à e (0, 1) e ((x% A) a sequência gerada pelo Algoritmo 
IALEM aplicado a EP (f, K), com o vetor erro associado e* e Rr. 

Seja ((x*, 9) a sequência gerada pelo Algoritmo IPPEMI aplicado a EP(L, R" x Rm), 
com vetor erro (e*, 0) e R" x R”, com os mesmos y, e ô. Se (xº, Aº) = (xº, A), então L(x*, A) 
= ((x*, N9), para todo k. 


Prova. A prova será feita por indução sobre k. Para k = O o resultado é válido por 
hipótese. Suponha que (x*, A) = (X%, N9. Precisamos mostrar que 


(1, NT) = (DM, NH). 
Em vista do Passo 2 do Algoritmo IPPEMI é necessário resolver EP (L, Rº x R7), 


com L£ e dada por (3.7), o qual tem solução única pela Proposição (3.5). Pela Proposição 
A 
(3.4) (es, A*) resolve o problema de minimização convexa definido como 
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Minimizar Fu à Se) (7, À) 


sa. (7,A) ER” x RZ. (3.14) 


com FO, A, A) =te (64, )9, (x, À)). As restrições deste problema são afim tal que 
a condição de Qualificação de Slater, dada na definição 3.2, se verifica e contudo existe um 
vetor KKT uí e R” (ver, por exemplo, Teorema 2.3.2 do Capítulo VII de (HIRIART- URRUTY; 
Lemaréchal, 1996), o qual trata com do caso não suave) tal que 


[ye(zt — 28) + ek] E [OF a (2E) + 5 MEhs(2*)], 


i=1 (3.15) 

NÓ) + vc AD=UL (I<ism), (3.16) 
A 

N>0, (3.17) 

u>0, (3.18) 

Mul=0 (1<i<m). (3.19) 


A A 
Observe que (3.15) e (3.16) são os subdiferenciais de L £ ((X8, A“, (x, À)) calculados 
A 
com respeito a xe À, respectivamente, no ponto (x, A)= (54, A”. Além disso, foi usado nestas 
equações, a hipótese de indução. Combinando (3.15) a (3.19) obtemos: 


hi(2*) 


XE = max (o X+ 
Yk 


h (1<i<m). (3.20) 
Substituindo (3.20) na equação (3.15) temos: 


[ye(m* — 28) + e] E [ofie(4) - = max (o A ur ohá)] E (3.21) 
k 


i=1 


No Passo 2 do Algoritmo IALEM, precisamos obter a solução X* de EP (is, R,). 
Dessa forma, aplicando a Proposição 3.2 obtemos que X* e EP Ea R,) se, e somente se, 


[me(zt — 28) + &] € [ofia(2) 1 q (o JE + dd ato) homo (z DE (3.22) 


i=1 


Usando a hipótese de indução, ou seja, x* = X* e A” = A4, obtemos de (3.21) que 
(3.22) se verifica com X* no lugar de $*e, então, X também resolve EP( (Le, R”). Desde que, 
pela Proposição 3.5, este problema tem solução única, podemos concluir que 
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X= +. (3.23) 


l 
Pelo Passo 4 do Algoritmo IALEM, tem-se x*'= 54 — — ek e, por outro lado, pelo 
1 do 
Passo y, do Algoritmo IPPEMI, temos x! = s* — ri o que pode ser concluído, pela 
k 
equação (3.23), que x! = x*!, Ainda falta a conclusão em relação à variável dual. Pelo 
Passo 3 do Algoritmo IPPEMI À“! = A4. Pela equação (3.20) 


ha(g* ) 
Yk 
(3.24) 


A+ = SE = max (o X + h (I<i<m). (3.24) 


Por outro lado, pelo Passo 3 do Algoritmo IALEM, tem se 


hi(5*) 


At = max (o + 
Yk 


h A<i<m). (3.25) 


Utilizando a hipótese de indução, A“ = A“, em (3.25) juntamente com (3.23) obtemos que 
(3.24) é igual a (3.25), ou seja A**! = A, completando o passo indutivo e a prova. 
A próxima proposição trata da terminação finita do Algoritmo IALEM. 


Proposição 3.6 Suponha que o Algoritmo IALEM pare na iteração k. Então, o vetor 


x“ gerado pelo algoritmo é uma solução de EP (f, K). 


Prova. Ver (NASRI, 2010) 
Finalmente, com o uso do Teorema 3.2 a prova de convergência do Algoritmo IALEM 


pode ser completada, pelo seguinte teorema: 
Teorema 3.3 Considere EP (f, K). Suponha que 


i) f satisfaça P1-P5 em Rº x Rº, 

ii) Ké dado por (3.1), 

iii) A Condição de Qualificação de Slater dada na definição 3.2 se verifica para K e 

iv) y, < (8, À] para algum À > 6, sendo 6 a constante de sub monotonicidade de f 
dada em P4. 


Considere f(x*, A”)) a sequência gerada pelo Algoritmo IALEM para resolver EP (f, 
K). Se EP (f, K) tem solução, então a sequência ((x*, N9)) converge para algum par ótimo 
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(x*, A*) para EP (f, K) e, consequentemente, x* e S(f, K). 
Prova. Ver (NASRI, 2010) 


É importante notar que a maior dificuldade a ser superada na implementação 
do Algoritmo IALEM é determinar a solução do subproblema gerado no Passo 2. 
Da demonstração do Teorema de Equivalência 3.2 ficou evidente que X* e S(L%., R) se, e 
somente se, a equação (3.22) é satisfeita, ou seja: 


m 


k ils k ha(ã) 
ev E (dFa(10) + > max 0, A; + —— 
i=1 


h on,(3%) + ya(E — ab), (3.26) 
Yk 


na qual dFX* (X*) é o subgradiente de F, calculado em X*, com F, dada por (3.9), oh, 


(X*) é o subgradiente de h, calculado em x*. Contudo, a Proposição 3.6 diz que se existe um 
procedimento que converge para a solução do problema de equilíbrio irrestrito EP (Ls, Rº) 
com e“ = 0. Então tal procedimento, após um número finito de passos, irá gerar um par (X*, 
e”) satisfazendo a condição do Passo 2 daquele algoritmo. 

Cabe agora uma explicação do real significado do vetor erro e* nos Algoritmos 
IALEM, IPPEM e IPPEMI, detalhada em (NASRI, 2010) e resumida a seguir. 

Suponha que temos uma sub-rotina que é usada para resolver o k-ésimo 
subproblema gerado no Passo 2 do Algoritmo IALEM com e* = O. Neste cenário, a sub- 
rotina gera soluções aproximadas X'*'cujos vetores erros são e*“' para /= 1,2, Uma vez 
que X*! é calculado em cada iteração / da sub-rotina, esta verifica se o vetor erro resultante 
e“! satisfaz o critério dado no Passo 2. Se assim for, o par (X*', e“) é aceito e o algoritmo 
inicia uma nova iteração k = k + 1. Caso contrário, o algoritmo executa a sub-rotina para 
mais uma iteração e define /= | + 1. Em particular, considere o caso quando fe h, são 
diferenciáveis. Nesta situação, usando (3.26), x *' satisfaz 


et = (5%! — 2) + VEge(56!) + s max (o xs da Vhi(55*). (3.27) 
3=1 ú 

Se x“! é a solução exata do k-ésimo subproblema, então o lado direito de (3.27) 
desaparece. Caso contrário, X*!' aproxima a solução exata e o lado direito de (3.27) é 
diferente de zero e é chamado de e*'. Agora, se e*' satisfaz a desigualdade do erro requerida 
no Passo 2 do Algoritmo IALEM com X*= x%!!, então consideramos e* = e*! e executamos 
o Passo 3 do algoritmo. Se x*' proposto não é bom o suficiente, é necessário um passo 
adicional da sub-rotina auxiliar, após o qual o teste será repetido com X* = xt"! e et= es, 
Pela Proposição (3.6), após um número finito de iterações, um par adequado (x'*, e”) deve 
ser gerado satisfazendo a desigualdade do erro requerida no Passo 2 deste algoritmo. 

No caso suave, cada iteração k do Algoritmo IALEM consiste em resolver uma 
equação não linear G(x) = 0, em que G: R" > R" é dada por 
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G*(a) = qa(x — a) + V Fila) + > max (o NM + 


i= 


fila) ! Vhia). (3.28) 


k 


Observamos que o lado direito de (3.28) é contínuo, porém não diferenciável, devido 
à presença do máximo. Em (FERREIRA, 2013) foram propostas duas metodologias para 
resolver o problema de equilíbrio irrestrito utilizando a relação (3.28) ao invés do Passo 2 
do Algoritmo IALEM. Uma das abordagens utilizadas substituiu o termo não diferenciável 
por uma aproximação suave que o tornou diferenciável e, então, o método de Newton foi 
utilizado para resolver o sistema suave gerado. A outra abordagem consistiu na utilização 
de uma família de métodos subgradientes para minimizar || G“(x)II? com G* dada por (3.28). 

Dando continuidade aos trabalhos desenvolvidos por (FERREIRA, 2013) e (NASRI, 
2010), esta pesquisa tem como um dos principais objetivos a resolução do problema de 
equilíbrio utilizando o método de lagrangeano aumentado, porém com penalidade diferente 
daquela usada em (NASRI, 2010). Desta forma, o subproblema gerado no Passo 2 do 
Algoritmo IALEM, será equivalente a um sistema não linear do tipo (3.28), e então, podemos 
utilizar métodos de otimização diferenciável para resolvê-lo. 
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LAGRANGEANO AUMENTADO EXPONENCIAL PARA O 
PROBLEMA DE EQUILÍBRIO 


Neste capítulo apresentaremos a principal contribuição desta tese, ou seja, o Método 
de Lagrangeano Aumentado Extragradiente Exponencial Inexato aplicado ao Problema 
de Equilíbrio EP (f, K) dado na Definição 3.1. O método de lagrangeano aumentado 
exponencial, aplicado ao problema de otimização, foi descrito no Capítulo 2 e agora será 
estendido para o Problema de Equilíbrio Geral, como desenvolvido em (NASRI, 2010) e 
(FERREIRA, 2013). Além disso, abordaremos aspectos teóricos e computacionais acerca 
do algoritmo proposto. 


11 DESCRIÇÃO DOS ALGORITMOS 


Nesta seção, apresentaremos o algoritmo de lagrangeano aumentado extragradiente 
exponencial inexato, assim como o algoritmo de ponto proximal extragradiente inexato. 

Considere p,: IR" x IR" x IR” x IR, — IR, parai=1,2,...,m,e a função lagrangeano 
aumentado L : IR" x IR" x IR”? x IR" x IR,.— IR dadas, respectivamente, por 


pi(2,4,A,7) = 1 [uertitO — xerti(] (4.1) 
e 
Hey) Hey)+ne-zy-z) +53 piz,y,A, 7). (4.2) 
i=1 


Neste caso, o Algoritmo IALEM 3.1 será reescrito incluindo estas novas funções e 
será chamado de IALEMe (o último e é de exponencial). 


Algoritmo 4.1 Lagrangeano Aumentado Extragradiente Exponencial Inexato 
(IALEMe) 


Passo 1: Considere uma sequência limitada ty ) c R,.. (0,1) e Rºx R”, uma 
tolerância de (0,1) ek=0. 


Passo 2: Encontrar um par (x*, e*) tal que X* resolve EP (8, R”) com ie definida 


como: 


Lea, y) = f(a, y) E Yale e E yY— LT) + Spitz, y, x, Yk) = (et, y— 2) 
i=1 


com p, dada por (4.1) e e* dado por: 
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les < dy IX = x“. 
Passo 3: Atualize a variável dual 


1 pi(ã : 
» ada = Me til , (= Last): 


Passo 4: Se (x*, A*) = (x *, A*!), então pare. Senão faça: 


k=k+1 
e volte ao Passo 2. 


Da mesma forma o Algoritmo IPPEMI 3.3 será reescrito e será denominado IPPEMle. 
NOTA: será usada a notação vetorial para representar o logaritmo natural de u e 
R” 


++ 


ou seja, In(u) = (In(u,), .. ., In(u,))”. 
Algoritmo 4.2 Ponto Proximal Extragradiente Inexato (IPPEMIe) 


Passo 1: Considere uma sequência limitadaly pc R,, (O, A)eR'xR”  ,eclR, 
uma tolerânciad e (0,1) ek=0. 


Passo 2: Encontre (X%, N) e Rrx R7  e(e% 0) e IR" x IR” tal que (X*, N9 resolve EP 
(Le, IR" x IR”) com 


Le (06, A), (y 4)) = Lx, A), (2) +y Axo, y —x)+y 4In(A)-In(A9, v-A)—(e%, y —0), 
(4.3) 
L dada por (3.2) e e* satisfazendo 


(e%, O) = les < dy (XE = x, Me N9I 


Passo 3: Se (es, AS = (x&, 9, então pare. Senão faça: 
l 


xt = k— e AS = & 
Vk 


k=k+1 evolte ao Passo 2. 


Observe que a diferença fundamental entre os Algoritmos IPPEMIe e IPPEMI 
está na função que regulariza o problema. Na função (4.3), o termo de regularização é 
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vAIn(A) —In(A9, 4 — À) ao invés de y, (A — A&, y — A). Com esta função é possível mostrar a 
convergência do Algoritmo IPPEMle, assim como a equivalência (dual) entre os Algoritmos 
IPPEMlIe e IALEMe. 

A seguir apresentaremos a convergência do Algoritmo IPPEMlIe. Estamos supondo 
que a boa definição do Algoritmo IPPEMIe é decorrente da boa definição do Algoritmo 
IPPEMI, tomando-se os devidos cuidados com o domínio da variável dual À em In(A), ou 
seja, À deve ser não negativo. 


21 CONVERGÊNCIA DO ALGORITMO PROPOSTO 


Nesta seção, nosso objetivo é tratar de aspectos teóricos relativos à convergência do 
Algoritmo 4.2, que estamos propondo. Conjecturamos que o mesmo é convergente. Porém, 
mostramos apenas uma espécie de convergência local, ou seja, conseguida impondo uma 
condição sobre o vetor de multiplicadores de Lagrange. 

O resultado seguinte é clássico considerando a literatura de análise convexa, ver 
por exemplo (SOLODOV, 2005). 


Teorema 4.1 Considere D c R” um conjunto aberto e convexo e ft: D >» R uma 
função diferenciável em D. Então as propriedades seguintes são equivalentes: 


a. Afunção fé convexa em D. 
b. Paratodosx,ye D, 

H(y) = f09) + (VF(X), y— X). 
c. Paratodosx, ye D, 


(VE) - VE), y- x) >0. 


Quando f é duas vezes diferenciável em D, as propriedades acima também são 
equivalentes à 
d. Amatriz Hessiana de fé semidefinida positiva em todo ponto de D, ou seja 


(Pra, d)>0, vxE D vde Rº. 


Prova. Ver (SOLODOV, 2005). 
Destacamos que no caso da função ser côncava as desigualdades acima continuam 
verdadeiras mas com sinais trocados. Por exemplo, no item (a) teremos, para todo x, ye D, 
É) <09 + (VÍ 09, y — 2). 


Lema 4.1 Considere f: R,, — R a função logaritmo natural f(x) = In(x). Então 
(i) fé uma função côncava. 
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(ii) In(x) <x— 1, Vxe Rj. 


Prova. A prova do item (i) é imediata. Basta derivar duas vezes a função fe usar O 
item (d) do Teorema 4.1. Para o item (ii), considere o item (b) do Teorema 4.1 juntamente 
com o item ()) deste Lema, ou seja, In(x) < In(y) + 8-0. Como a desigualdade vale para 
todo x, ye R,, basta escolher y = 1 e ter-se-ã completado a prova. 


Lema 4.2 Considereue R, ev we R. Então 
(1 +In(u)-u(w-v)z2065 w=<v. 
Prova. Desde que In(u) é uma função côncava, pelo item (ii) do Lema 4.1 
In(u)su-1teS1+In(u)-u<o. 
Consequentemente, 
1d +In(u-u(w-v)z06w-v<0S w<vyvyweRr. 
Teorema 4.2 Considereu, ve R”,, eweR”. Então, 


(v-u-In(y+Inu),w-)z06w=<sv. 


Prova. (v- u-In(v) + In(u), w-— v) = Em [(v,— u) + (In(u) — In(vDI(w, — v) = Zm A 
i=1 i= 


(1 + In(u) — u)(w, — v). Na última desigualdade foi usado o item (ii) do Lema 4.1, ou seja 
In(v)sv-1S-In(v)+vz1,vVi=1,...,m. Pelo Lema 4.2, para cada i=1,...,m, (1 + 
In(u) —- u)(w- v)z0S w<v, vw e Re portanto a prova está completa. 

A seguir provaremos a convergência do Algoritmo IPPEMlIe, no sentido do que foi 
explicado no início da seção, prova esta que será baseada na convergência do Algoritmo 


IPPEMI. 


Teorema 4.3 Convergência do Algoritmo IPPEMIe 

Considere EP (f, K) com K dado por (3.1), f satisfazendo P1-P5 em R" x R" e uma 
sequência ly) < (0,Y] para algum y > 6, em que 6 é a constante de sub monotonicidade 
de f. Seja “x k À 9 a sequência gerada pelo Algoritmo IPPEMIe com À“ > u>0. Seo 
problema EP (L, R" x R”,,) tem solução, então 0%, ÀS) converge para algum par (x*, A*) 
esS(L, RºxR7). 
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Prova. Considere as sequências xs, N 9jc Rºrxa” e(e%c R” geradas pelo 
Algoritmo IPPEMle, 4.2, então pelo Passo 2 tem-se 


Le((8, À, (7, ))>0 (4.4) 
e 

Ile < dy UIÇÃ = 8, MAI (4.5) 
em que 
L(06,Ã9, (1,40) = UQO, AY, (4 40) + VU — 8, y — + (4.6) 


+YAIn(ÃO — In(A9, 4 — À — (65 y— 5, 


e L é dada por (3.2). 
Além disso, pelo Passo 3 do Algoritmo IPPEMlIe tem-se 


1 ' 
pro. o Nº=r (4.7) 


Por outro lado, substituindo estes valores na função Le e também na fórmula do erro 
do Passo 2 do Algoritmo IPPEMI, 3.3, tem-se 


LS((2, A), (y, 10) = L((28, E), (y, 10) + (2t at, yo) Ap (DEAR, DB) — (ey 2), 
(4.8) 


com L dada por (3.2). Já a fórmula do erro coincide com (4.5) e, da mesma forma, a 
atualização de x e À no Passo 3 de IALEMI coincide com (4.7). Agora, subtraindo (4.6) de 
(4.8) temos: 

a A A A A A A A 
LE (6 AY, (4 2)) — E (0E AY, (0440) = y4ln(ÃO = In(A, e AM — y4N E = AS, o AM. 
Isolando Leno primeiro membro desta última igualdade e usando (4.4) tem-se 

Le (GOÃO, (410) 2 y 8 AS = A 9 = y4n(ÃO — In(A9, pr — ÀS) 
ou, equivalentemente 
Le (08, Ã9, (4,40) 2 y AS AS= In(Ã£) + In(A9, = ÀS. (4.9) 
Por hipótese yu < Às então pelo Teorema 4.2 tem-se que 


(e 8 = In(Ã9 + In(A9, /— À9 >0. 


Desde que y, >0, então (4.9) se torna. 
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L(0%, 9, (1) 20 (4.10) 
Logo, de (4.10), (4.5) e (4.7) tem-se todos os passos do Algoritmo IPPEMI satisfeitos 
e, portanto, pelo Corolário 3.1 se EP (L, Rº x R” ) tem solução então “xs, À 9) converge 
para algum par (x*, A“) E S(L, R" x R”,)e a prova está completa. 


31 EQUIVALÊNCIA ENTRE IALEME E IPPEMIE 


Nesta seção, mostraremos a equivalência entre as sequências geradas pelos 
Algoritmos IALEMe e IPPEMlIe e, desta forma, obtemos a convergência do Algoritmo 
IALEMe. 


Teorema 4.4 Equivalência entre IALEMe e IPPEMIe Considere EP (f, K) satisfazendo 
P1 - P4 em Rº xRº, uma sequência (y,) e R,, um erro relativo de tolerância d e (0, 1) e 
((x*, AM)y a sequência gerada pelo Algoritmo IALEMe aplicado a EP (f, K), com vetor erro 
associado e* e R". Por outro lado, considere ((x*, A)y a sequência gerada pelo Algoritmo 
IPPEMle aplicado ao EP(L, R"x RR”), comL dada por (3.2) e vetor erro associado (e*, 0) e 
R" x Rr, usando os mesmos y* e 5. Se (xº,Aº) = (Xº, x), então ((x*, A9)y = L(x*, A) para todo k. 


Prova. A prova é feita por indução sobre k. Para k = 0, o resultado é verdadeiro por 
hipótese. Suponha que (x*, A”) = (X*, A”) (Hipótese de indução). Provaremos que 


(xt As?) mi (e Ne). 
Em vista do Passo 2 do Algoritmo IPPEMlIe é necessário resolver EP (L£, R" x R”) 
A 
com L2:RºxR, — R dada por (4.3). Pela Proposição (3.4), s(6S, A) resolve o problema 


de minimização convexa definido como 


dao O AA A 
Minimizar F(X*, N*)(x, À) 


sa (x Ae Rºx Rm. (4.11) 
A 
com FS, IN; x, À) = Ls ((X%, À“), (x, A), sendo Ls dada por (4.3). Observe que, 


AA A A, A A 
o,F(xS A) (x, A) = Le ((x8 AM, (x, À)) = OL (XS AS), (x, À) + y txt — x) — es, 


com L(($5 M), (x, A) = (4%, x) + 3 Mni(ax) — 3 Mhi(2º) (ver relação (3.2)) e 
2=1 1=1 


IL (68, NE) (x, 1) = df(s + DNtOhi(a). 


i=1 
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No ponto x = ** e usando F&y) = f(x, y) esta última equação resulta 


Or L((2*, AP), (28,2) = 0Fa (28) + 3 XfOhi("). Logo, 0 E OF (se su)(28,2) & 
i=1 


0€ [nst-at)-é + 0Fos(£ 43X Ohil en) e 


1=1 
é [yu(at — 2%) + e] € (oro + Dito 2) 
Estas equações serão utilizadas nas condições KKT para o problema (4.11). Desde 
que as restrições do problema (4.11) sejam afins, a Condição de Qualificação de Slater 
dada na Definição 3.2 se verifica. Contudo existe um vetor de KKT v“e R” (ver, por exemplo, 


Teorema 2.3.2 do Capítulo VII de (HIRIART-URRUTY; Lemaréchal, 1996), o qual trata com 
o caso não suave), tal que 


[elx — 2º) + e] E (oria(i) +43  Mohi(ê 5) (4.12) 
il 
NB +y ANAC InÃD=Vs  di=1,..,m (4.13) 
A 
Mt20, d=1,...,m (414) 
V20, d=1,...,m (4.15) 
A 
My=0,i=1,...,m (4.16) 


Observe que (4.12) e (4.13) são os subdiferenciais de Le (8, NJ (x, À)) calculados 
com respeito a xe À, respectivamente, no ponto (x, A)= (e, NS Além disso, usamos nestas 
equações, a hipótese de indução. 

De (4.13) - (4.16) temos 


A = AZ 1a 
me [ln dE — Im AB] = hi(28) 6 In(Ã8/28) = le!) 
k 
ou seja, 
dE — Abe tét), (4.17) 


Substituindo (4.17) em (4.12), obtemos 


Into at) +] € [ora(i! 45 MEME an (88), i=1,.. om (418) 
1=1 
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No Passo 2 do Algoritmo IALEMe precisamos determinar a solução X* de EP (Ls, Rº) 
com Ee R”x Rº > R dada por 


Le(x,y) = f(x,y) + 9yulz — q*,y— z) + > [ee (9) — Ate? e ha eo) (&,y — 2). 
Para isso, aplicando a Proposição 3.2 a este problema temos que X* e s(Ls , Rº) 
se, e somente se, 


Jehi(Et) 


[rue (a — 8) + e E [9fia(a) + = Men ons(8) : (4.19) 


Pela hipótese de indução x = x*e X*= À 4 então (4.19) torna-se 


[rue (7* — q +] e ora! 45 He a 


i=1 


onda! . (4.20) 


Desde de que, pela Proposição 3.5 este problema tem solução única, então (4.20) 
coincide com (4.18) quando 


Xt= x (4.21) 
Pela atualização de x, no Passo 4 do Algoritmo IALEMe, 


= =p (4.22) 


1 
Ft o -&. (4.23) 


De (4.21), (4.22) e (4.23), concluimos que x! = x 1, 
Pela atualização das variáveis duais no Passo 3 do Algoritmo IPPEMIe e por (4.17), 
temos: 


E+ = Moment, (4.24) 


se : 
Agora, comparando (4.24) com At = a : dado no Passo 3 do Algoritmo 
4.1, juntamente com a hipótese de indução, concluímos que A**! = À *!, completando o 


passo indutivo e a prova. 
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41 AJUSTE QUADRÁTICO DA PENALIDADE EXPONENCIAL 

Na seção 3.2, em que foi considerado o caso suave, a cada iteração k do Algoritmo 
IALEM era necessário a resolução de um sistema não linear G(x) =0, em que G:R" > Rº 
é dada por 


G*(x) = (x — 2º) + V Fi(x) 4 mar (o XE + aa Vhi(x). (4.25) 
i=1 k 


A principal diferença entre os Algoritmos IALEM e IALEMe é uso da penalidade 
exponencial no lugar da penalidade quadrática clássica. Desta forma, com relação ao 
Algoritmo IALEMe, o sistema não linear 4.25 pode ser reescrito como 


GH(2) = (x — 2º) + VFi(2) +45) Me yh(a). (4.26) 


i=l 

Em geral, os métodos de lagrangeanos aumentados exponenciais tendem a gerar 
subproblemas mau condicionados devido ao rápido crescimento da função penalidade 
exponencial. Uma alternativa para contornar o mau condicionamento é fazer um ajuste 
quadrático da função penalidade, como utilizado em (MARTINEZ, 2000; MATIOLI, 2001), 
entre outros. 

Adotando esta estratégia faremos um ajuste quadrático da função exponencial. 
Considere À e R,, fixado e a função p: R — R dada por p(y) = Ae ep:R-> R dada por 


y 
Ae eye 
-| ai is (4.27) 


Pg) = ay +by+c sey>o0. 
O objetivo do ajuste é garantir a suavidade da função, da primeira e segunda 


derivadas de p. Denotando g(y) = ay? + by + ce p(y) = Ae ; às partes quadrática e 
exponencial, respectivamente, de p, serão impostas as condições de interpolação: 


p(O) = q(0), p'(0) = q(0) e p4(0) = q(0). 


Da condição de interpolação p(0) = g(0) temos À = c. A derivada da função p é dada 
por 


AA <0 
E Dem  sey< 
P(y)=4 % (4.28) 


2ay+b sey>o), 


e usando a condição de interpolação da derivada de p, p'(0) = g'(0), temos que b 


Bl» 


A segunda derivada da função p é dada por 
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À yu 
=” Sem sey<0 
Pl) =4 % (4.29) 


2a se y > 0. 


Assim, usando a condição de interpolação na segunda derivada de p, p'“(0) = q“(0) 


temos que a = —s. 
2 
2Y%k 


E 
Logo, q(y) = 23! + da A e consequentemente, 
k 
E» 
Ae7k sey<0 
(4.30) 


A A 
3 +—y+A sey>0 
2% Ye 


Com a quadrática ajustada, o sistema não linear (4.26) torna-se 
Ga) = (x — 2º) + V Fi(2) +55 Dilh(2))V hi(x). (4.31) 
i=1 


em que p, parai=1,...,m, é dada por (4.30). 
No próximo capítulo faremos testes numéricos envolvendo o algoritmo proposto, 


IALEMe, com e sem o ajuste quadrático da penalidade exponencial. 
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EXPERIMENTOS NUMÉRICOS 


11 INTRODUÇÃO 


Com o intuito de validar a nova metodologia, apresentaremos a seguir resultados 
numéricos obtidos com a implementação do algoritmo proposto. 

O Algoritmo IALEMe, 4.1, foi implementado em MatLab versão 8.1.0.604 (R2013a), 
executado em um computador MacBook Air com processador 1.4 Intel Core i5 e utilizado na 
resolução de problemas de equilíbrio geral e problemas de equilíbrio de Nash Generalizado 
(GNPs). A classe de problemas desiganados GNEPs, sigla derivada de sua denominação 
em inglês, compõe uma classe particular do problema de equilíbrio geral. 

Para iniciar, escreveremos o problema de equilíbrio de Nash Generalizado em termos 
de um problema de equilíbrio geral. Para tal, descreveremos brevemente o GNEP. Para o 
estado da arte e uma descrição completa do problema de equilíbrio de Nash Generali- zado 
pode ser consultado, por exemplo, (FACCHINEI; KANZOW, 2010a; FACCHINEI; KANZOW, 
2007) e suas referências. 

O termo GNEP designa um jogo não cooperativo com N jogadores no qual cada 
jogador j controla sua variável de decisão, designada por x' e RY”. A variável de decisão de 
cada um dos outros jogadores é dada por 


O vetor x = (x, x) = (X,...,X1, x, X,...,xN) e X denota a variável de decisão 
N 
de todos os jogadores, em que X R" é um conjunto convexo fechado e n= > n; 
j=1 
Formalmente, resolver um GNEP equivale a resolver o problema de 


min 9 (x, x?) sa: xeX(x), (5.1) 
>»d 


onde X(x7) « R” e 6, denotam, respectivamente, o conjunto de estratégias e a 
função de retorno do jogador j para cada j= 1,..., MN. 
Apresentaremos, a seguir, exemplos nos quais o conjunto X tem a forma 


X=(xeR":w(xy)<O, Vi=1,...,G, (5.2) 


no qual w,: R” — R é uma função convexa para cada i= 1,2,...,4. 

Para a resolução destes problemas, inicialmente reformulamos o GNEP em termos 
de um problema de desigualdade variacional, como descrito em (FACCHINEI ANDREAS 
FISCHER, 2007) e, então, usando o fato que o problema de desigualdade variacional pode 
ser reformulado como um problema de equilíbrio geral (ver, por exemplo, (IUSEM; SOSA, 
2003a)), obtemos um problema de equilíbrio representado por um GNEP cuja função 
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resultante f: R"x R"—> Re o conjunto K—> R" têm, respectivamente, a forma 
N . . 
fle)=5 (Voa), y 2) eK=X. (5.3) 


Diante do exposto, o problema da poluição de uma bacia hidrográfica, a ser tratado 
na próxima seção, é formulado em termos de um GNEP. Contudo, seguindo a discussão 
acima, este problema pode ser tratado como um problema de equilíbrio. Vale destacar que 
este é um problema importante e que tem sido objeto de estudos resultando na publicação 
de diversos artigos. Entre estes, (KRAWCZYK; URYASEV, 2000; MATIOLI; SOSA; YUAN, 
2012; SANTOS; SCHEIMBERG, 2011), merecem destaque. 


21 APLICAÇÃO DO ALGORITMO IALEME COM PENALIDADE EXPONENCIAL 

Nesta seção apresentaremos exemplos onde aplicamos o algoritmo proposto na 
resolução de problemas de equilíbrio de Nash Generalizado e problemas de equilíbrio 
gerais. 

O problema da poluição de uma bacia hidrográfica, com 3 jogadores, e sua generaliza- 
ção para N jogadores é o primeiro a ser tratado. O outro, um problema de equilíbrio geral 
com 5 jogadores e sua generalização para N jogadores. Uma lista completa com estes e 
outros problemas pode ser encontrada, por exemplo, em (FACCHINEI; KANZOW, 2009). 


Exemplo 5.1 No problema da poluição de uma bacia hidrográfica com três jogadores 
j=1,2,3, estes figuram como possíveis agentes poluidores. O problema consiste em 
minimizar a função objetivo dada por 


9,()=c x2+d x(x + X2 +) - Dx, j=1,2;9, (5.4) 


emquec, =c,=0,01,c,=0,05,d,=0,01,b,=2,90,b,=2,88 andb,=2, 85. 
As funções w, dadas em (5.2), são definidas por 


w(M)=3,25x +1,25X+4,125x — 100 
wx) = 2, 291x! + 1, 5625Xº + 2, 8125xº — 100. 


Para construirmos o sistema não linear G4x) = 0, dada em (4.26), será necessário 
calcular VF (x), o qual é dado por 
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VFs(z) = | 0,017! 40,122? + 0,017º — 2,88 
[ 0,017! +0,0122 + 0,042º — 2,85 


| 0,042! + 0,0172 + 0,01x3 — 2,90 | 
| 


Neste exemplo, os valores iniciais escolhidos foram: X = [5,9,3]",Nº=[1,1.5]" e 
Y,= 0.4 para todo k. 

Na Tabela 1 a seguir, a primeira coluna representa o nome dos algoritmos utilizados 
na comparação. A segunda, contém a solução obtida pelos algoritmos comparados e 
a terceira e quarta colunas são, respectivamente, o número de iterações e o tempo 
computacional, em segundos, gasto na resolução do problema. 

Os algoritmos comparados são: 


* - IALEMe: Algoritmo de Lagrangeano Aumentado Extragradiente Exponencial 
Inexato. 


* RA:Algoritmo de Relaxação, descrito em (KRAWCZYK; URYASEV, 2000). 


*-— MSY: Um Algoritmo Numérico para Encontrar Soluções de um Problema de 
Equilíbrio de Nash Generalizado, descrito em (MATIOLI; SOSA; YUAN, 2012). 


* | SS: Um Algoritmo de Subgradiente Inexato para Problemas de Equilíbrio, des- 
crito em (SANTOS; SCHEIMBERG, 2011). 


*-— EUDA: Algoritmos para o Problema de Equilíbrio Aplicados ao Problema de 
Equilíbrio de Nash, descritos em (FERREIRA, 2013). 


Algoritmo Solução obtida lterações tempo 
IALEMe [21, 144; 16, 028;2, 7268] 125 0,09 
RA [21, 149; 16, 028; 2, 722] 20 E 
MSY [20, 115;17,125;8, 205] 10 1,37 
ss [21, 144; 16.028 2.726] 7 - 

EUDA [21,143;16, 028;2, 727]T 6 - 


TABELA 1 — Resultados do Exemplo 5.1 


Analisando os resultados da Tabela 1, podemos tirar algumas conclusões. Primeira- 
mente destacamos que algoritmo proposto, IALEMe, utilizou mais iterações para resolver 
o problema quando comparado com os demais algoritmos. Porém, dispendeu tempo 
computa- cional menor quando comparado com o algoritmo MSY. O motivo principal do 
Algoritmo IALEMe utilizar um número maior de iterações tem a ver com o fato da penalidade 
exponencial ser muito sensível a pequenas alterações nos dados. Isso se deve à presença 
do parâmetro y na expressão do sistema não linear, (4.26), a ser resolvido a cada passo do 


algoritmo, gerando, em geral, problemas mal condicionados. 
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Na sequência executamos o Algortimo IALEMe, aplicando-o ao mesmo problema, 
porém agora com quadrática ajustada. Os resultados são mostrados na Tabela 2 a seguir, 
com as colunas sendo as mesmas já descritas na Tabela 1. 


Algoritmo Solução obtida lterações | tempo 
IALEMe [21, 144;16, 028; 2, 7260] 3 0,04 
RA [21, 149;16, 028;2, 722] 20 - 
MSY [20, 115; 17,125;8, 205] 10 1,37 
ss [21, 144; 16, 028;2, 726] 7 E 

EUDA [21, 143; 16, 028;2, 727] 7 6 = 


TABELA 2 — Resultados do Exemplo 5.1 com quadrática ajustada 


Observe agora, nos resultados da Tabela 2, que houve grande melhoria no 
desempenho do Algoritmo IALEMe. 

O ajuste quadrático consegue contornar o problema do mau condicionamento, ou 
seja, o algoritmo ficou menos dependente do valor de y na expressão do sistema não linear 
em (4.26). Ainda, podemos observar que, o algoritmo apresenta desempenho melhor em 
relação aos demais com os quais foi comparado. 


Exemplo 5.2 Neste exemplo, consideramos a generalização do problema da poluição 
de uma bacia hidrográfica para N jogadores ao invés de 3. Os testes são executados 
variando-se o valor de N. Neste caso, a função 0 de cada jogador é dada por 


N 
8(x) = u(x)? + px! (e) VM Vj=1,..0N, 
3=1 


onde 8 = 0,01 eu, v,€ R para cada j=1,..., N. Com esta generalização temos 
m restrições, as quais representam as constantes de emissão de poluição, definidas como 
Ax<b, ondeAe R”Nebe R”. O vetor de multiplicadores de Lagrange Àº, a matriz A, o 
vetor b e os escalares u ev, são gerados randomicamente. O ponto inicial foi escolhido 
como sendo X =[1,1,...,1]".A função f associada ao problema de equilíbrio nesta 
configuração é dada por (5.3) e o conjunto Ké dado por K= X=(xe R": Ax< b). 

Na Tabela 3 a seguir, a primeira coluna representa o nome do algoritmo, a segunda 
o número de jogadores, a terceira o número de restrições, a quarta o número de iterações 
executadas e a última coluna representa o tempo, em segundos, gasto para resolver o 
problema em questão. 

Consiatamos que, para este problema, as respostas são melhores quando o 
parâmetro y é pequeno. Por isso, nos testes consideramos y, = 1e — 5. Além disso, para 
problemas maiores (acima de 1000 jogadores), o método não consegue obter uma solução 
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para o problema devido ao mau condicionamento gerado pela penalidade exponencial. 
Esta dificul- dade é contornada com a utilização da penalidade exponencial com quadrática 
ajustada. Apresentaremos os resultados desse procedimento nos testes da próxima seção. 


Algoritmo N m | lterações | Tempo 
IALEMe 10 5 1 0,03 
IALEMe 30 | 20 2 0,03 
IALEMe 50 | 10 2 0,04 
IALEMe 100 | 10 2 0,4 
IALEMe 200 | 20 2 0,12 


TABELA 3 — Resultados do Exemplo 5.2 


Como os dados desse exemplo foram gerados randomicamente não seria conveniente 
realizar comparações. De qualquer forma, este mesmo exemplo foi apresentado em (FER- 
REIRA, 2013; MATIOLI; SOSA; YUAN, 2012) e os resultados são parecidos. Mas ainda 
persiste o problema do mau condicionamento gerado pela penalidade exponencial e só 
conseguimos resolver problemas com dimensões não muito altas, com até 200 jogadores e 
20 restrições, como mostrado na tabela acima. Por isso, esse mesmo exemplo foi executado 
usando a penalidade com quadrática ajustada e com os mesmos dados iniciais, ou seja, 
geramos randomicamente os dados e salvamos em disco para executar e comparar as 
duas abordagens, com e sem quadrática ajustada. Os resultados são mostrados na Tabela 
4 a seguir. 


Algoritmo N m lterações | Tempo 
IALEMe 10 5 1 0,03 
IALEMe 30 20 2 0,03 
IALEMe 50 10 3 0,04 
IALEMe 100 10 4 0,06 
IALEMe 200 20 2 0,11 
IALEMe 1000 20 7 70,49 
IALEMe 1000 | 100 9 289,34 


TABELA 4 — Resultados do Exemplo 5.2 com quadrática ajustada 


Observe que, com o ajuste quadrático, foi possível resolver problemas maiores com 
até 1000 jogadores e 100 restrições, o que não ocorria sem o ajuste quadrático. 


Exemplo 5.3 Neste exemplo, consideramos um problema de equilíbrio quadrático 
EP(f, K), para i=1, 2, introduzido em em (SANTOS; SCHEIMBERG, 2011) e também em 
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(TRAN; DUNG; NGUYEN, 2008), onde 


5 
K=DeRr: > We -6 sk D,/8 1,000, 5) 
j=1 
e 
fx, )=(Px+Qy+qy-x), i=1,2, (5.5) 
na qual 
[e 2 00 il [Re 2 0 0 q 
[2360 00] [2360 00] 
=| 0 0 45 20), B=|0 0/45 2 0); 
|0 0 2 330] |O D 2330] 
[0 0 0 0 2] [o 0 0 0 3] 


[156 10 o | [1] 
[1 1,6 0 | | —2 | 
Q=| 0 0 1,5 ol eg=|-1. 
|O 1,5 o | | 9] 
Vo 0 0 E [aii 


Os valores iniciais considerados para este exemplo foram: 
X=[1,3,1,1,2]T, A =[1,1,1,1,1,1,1,1,1,1,1] ey=le-s. 


Observe que a restrição -5 <x'<5 é equivalenteax- 5 <0e-x- 5 <0, para todo 
j=1,...,5. Os resultados numéricos são mostrados nas Tabelas 5 and 6 a seguir. Os 
algoritmos comparados são: 


* - IALEMe: Algoritmo de Lagrangeano Aumentado Extragradiente Exponencial 
Inexato. 


* | EA: Algoritmo Extragradiente Estendido para o Problema de Equilíbrio, descrito 
em (TRAN; DUNG; NGUYEN, 2008). 


* | SS: Um Algoritmo de Subgradiente Inexato para Problemas de Equilíbrio, des- 
crito em (SANTOS; SCHEIMBERG, 2011). 


Algoritmo Solução obtida lterações Tempo 
IALEMe [-0,725;0,803;0,720;-0, 867;0, 250] 1 0,04 
EA-a [-0, 726; 0, 804;0,720;-0, 866;0, 250] 10 - 

ss - 10 0,0006 


TABELA 5 — Resultados do Exemplo 5.3 associado com EP (f,; K) 
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Algoritmo Solução obtida lterações | Tempo 
IALEMe [-0, 725;0,803;0,720;-0, 867;0, 200] 1 0,04 
EA-b [-0, 726; 0, 803;0,718;-0, 864;0, 200] 21 - 

ss - 10 0,0006 


TABELA 6 — Resultados do Exemplo 5.3 associado com EP (f,; K) 


Para este exemplo, os autores (TRAN; DUNG; NGUYEN, 2008) não exibem o tempo 
computacional para o Algoritmo EA, para nenhum dos dois casos a e b. Se compararmos 
somente o número de iterações, o Algoritmo IALEMe foi superior. Por outro lado, os autores 
(SANTOS; SCHEIMBERG, 2011) não exibem a solução para o Algoritmo SS. Comparando 
então, o número de iterações, o Algoritmo IALEMe foi superior ao Algoritmo SS, porém o 


inverso ocorreu em relação ao tempo computacional. 


Exemplo 5.4 Neste exemplo, consideramos a generalização do problema anterior 
para N jogadores ao invés de 5. Desta forma, o conjunto K é dado por 
N 
K=(xeRN:S xz-1,-5<x 5, je 1t,..., N) 
j=1 
O número de restrições ém = 2N + 1, as matrizes P, parai=1,2, Qe o vetor q são 
gerados randomicamente. Os resultados numéricos são mostrados na Table 7 a seguir, 
com o parâmetro y, = te — 5 para toda iteração k. 


Algoritmo N m lterações | Tempo 
IALEMe 10 21 1 0,04 
IALEMe 20 41 1 0,04 
IALEMe 50 101 1 0,06 
IALEMe 100 | 201 1 0,09 
IALEMe 150 301 1 0,17 


TABELA 7 — Resultados do Exemplo 5.4 


Aqui também, os dados foram gerados randomicamente, e por isso as comparações 
com outros métodos não seriam convenientes. Mas, como ocorreu nos exemplos anteriores, 
quando utilizamos a penalidade exponencial sem o ajuste quadrático, só conseguimos 
resolver o problema com até 200 jogadores e 20 restrições. 

Na próxima tabela, listamos os resultados alcançados pelo Algoritmo IALEMe, mas 
agora com quadrática ajustada. Neste caso, para contrastarmos as duas abordagens, 
manti- vemos o valor de y, = te — 5 para os testes com até 200 jogadores e 20 restrições, 
que foram aqueles que o Algoritmo IALEMe, sem o ajuste quadrático, conseguiu resolver. 
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Para os demais testes, ou seja, acima de 200 jogadores, o Algoritmo IALEMe se comportou 
melhor quando y, é maior. Por isso adotamos o valor y, = 5, para todas as iterações k 
(apesar de resolver também com outros valores para y, inclusive o valor y, = 1e — 5, para 
todas as iterações k mas com desempenho um pouco pior). 


Algoritmo N m lterações Tempo 
IALEMe 10 21 1 0,04 
IALEMe 20 41 1 0,04 
IALEMe 50 101 1 0,06 
IALEMe 100 201 1 0,08 
IALEMe 150 301 1 0,09 
IALEMe 300 601 3 0,91 
IALEMe 500 1001 3 1,5 
IALEMe 1000 2001 3 11,2 
IALEMe 2000 | 4001 3 83,8 
IALEMe 3000 | 6001 3 291 
IALEMe 4000 | 8001 3 683 


TABELA 8 — Resultados do Exemplo 5.4 com quadrática ajustada 
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CONSIDERAÇÕES FINAIS 


Inicialmente procedemos a uma revisão bibliográfica dos tópicos de otimização 
necessários para o desenvolvimento da proposta desta tese. Destacamos como conceitos 
principais aqueles relacionados aos métodos de lagrangeano aumentado e ao problema de 
equilíbrio geral. 

Em um primeiro momento pretendíamos avançar no sentido de estender os algoritmos 
propostos em (NASRI, 2010) e (FERREIRA, 2013) para resolver o problema de equilíbrio 
geral. Cumprimos este objetivo propondo um novo algoritmo de lagrangeano aumentado 
com penalidade exponencial, em substituição à penalidade clássica de Rockafellar que 
estes autores utilizaram em seus trabalhos. Com esta abordagem introduzimos o Algoritmo 
IALEMe, o qual sintetizou a principal contribuição deste trabalho. 

Na sequência, procuramos refazer a teoria geral em torno do novo algoritmo. Mostra- 
mos que a sequência de multiplicadores (dual), gerada pelo Algoritmo IALEMe, coincide 
com a sequência gerada pelo Algoritmo IPPEMlIe com distância de Bregman, assim como 
foi realizado em (NASRI, 2010). 

Sobre a convergência conjecturamos que o algoritmo é convergente. Conseguimos 
mostrar uma espécie de convergência local impondo uma condição sobre o vetor de 
multiplicadores de Lagrange. 

Em um segundo momento, com o objetivo de contornar os problemas de mau 
condicionamento decorrentes do uso da penalidade exponencial, realizamos um ajuste 
quadrático desta. Realizamos testes numéricos empregando a penalidade exponencial e a 
também utilizando esta com quadrática ajustada. 

O algoritmo proposto foi implementado em MatLab e aplicado na resolução de 
problemas de equilíbrio e problemas de equilíbrio de Nash Generalizados, já conhecidos 
na literatura. Comparamos os resultados alcançados pelo IALEMe com outros trabalhos já 
publicados. 

Os testes numéricos mostraram que, utilizando o ajuste quadrático na penalidade 
exponencial, o IALEMe teve desempenho superior aos algoritmos com os quais foi 
compa- rado. Por outro lado, quando utilizamos a penalidade exponencial sem o ajuste 
quadrático, o algoritmo pode gerar subproblemas mal condicionados, ou seja, dependente 
do parâmetro de penalidade. No entanto, a penalidade exponencial foi importante para 
mostrar os principais resultados da tese. 

Como sugestões para trabalhos futuros, destacamos: 


* Apresentar novos algoritmos com penalidades diferentes da exponencial, tal 
como a logarítmica quadrática, com a qual tem-se obtido algum sucesso em 
trabalhos envolvendo o método de lagrangeano aumenado para o problema de 
otimização. Acreditamos que o mesmo possa ocorrer com o problema de equi- 
líbrio geral, desde que o problema de otimização pode ser reformulado como 
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um caso particular deste. 


Realizar testes numéricos com outros problemas de equilíbrio geral, incluindo 
casos particulares destes, como por exemplo, desigualdade variacional. 


Pesquisar estratégias diferentes para a escolha do parâmetro de penalidade y 
tais como aquelas desenvolvidas em (BIRGIN; MARTINEZ, 2012). 
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